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ABSTRACT 


Modern high accuracy measurements of the non-rigid earth are to 
be referred to four-dimensional, i.e., time- and space-dependent, 
reference frames. Geodynamic phenomena derived from these measure- 
ments are to be desc-ibed in a terrestrial reference frame in which both 
space- and time-like variations can be monitored. Existing conven- 
tional terrestrial reference frames (e.g. CIO, BIH) are no longer suit- 
able for such purposes. 

The ultimate goal of this study is the establishment of a 
reference frame, moving with the earth in some average sense, in which 
the geometric and dynamic behavior of the earth can be monitored, and 
whose motion with respect to inertial space can also be determined. 

The study is conducted in three parts. In the first part prob- 
lems related to reference directions are investigated, the second part 
deals with the reference origins and the third part with problems 
related to scale. 

The approach is based on the fact that reference directions at 
an observation point on the earth surface are defined by fundamental 
vectors (gravity, earth rotation, etc.), both space and time variant. 
These reference directions are interrelated by angular parameters, also 
derived from the fundamental vectors. The interrelationships between 
these space- and time-variant angular parameters are illustrated in 
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hierarchic structures or towers, which make the derivations of the 
various relationships convenient. In order to determine the above 
parameters from observations using least squares techniques, model 
towers of triads are also presented to allow the formation of linear 
observation equations. Although the model towers are also space and 
time variant, their variations are described by adopted parameters 
representing our current knowledge of the earth. 

After the translational and rotational degrees of freedom 
(origin and orientation) have been discussed, the notion of a length, 
scale degrees of freedom are introduced and studied under spacelike/ 
timelike variations. 

According to the notion of scale parallelism, originated by 
H. Weyl, scale factors with respect to a unit length are given. 
Three-dimensional geodesy is constructed from the set of three base 
vectors (gravity, earth- rotation and the ecliptic normal vector). 

Space and time variations are given with respect to a polar and singular 
value decomposition or in terms of changes in translation, rotation, 
deformation (shear, dilatation or angular and scale distortions). 
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PART I: SYSTEM OF REFERENCE DIRECTIONS : THE E-TOWER 


by 

Erik W. Grafarend, Ivan I. Mueller, 
Haim B. Papo and Burghard Richter 


1 


Introduction 


In order to take full advantage of high quality geodetic obser- 
vational systems, such as lunar and satellite laser ranging and radio 
interferometry to quasars, an appropriate terrestrial reference frame 
is needed in which geodynamic phenomena can be detected and monitoied. 
The importance of the definition, determination and subsequent mainte- 
nance of such a terrestrial reference frame has been recognized by 
many, although, so far, no satisfactory and comprehensive proposals for 
its realization have been put forward [Kolaczek and Weiffenbach, 1975; 
IAU, in press j. 

The ultimate goal of this study is the establishment of such a 
reference frame, moving with the earth in some average sense, and whose 
motion with respect to inertial space can also be determined. 

In attempting a solution to the problem, a "zero base" approach 
is taken. Being fully aware of the large body of accumulated knowledge 
in the relevant disciplines of geodesy, astronomy and geophysics, we 
conduct a step-by-step analysis of known concepts and relationships with 
the purpose of establishing an unbiased and systematic foundation. In 
many cases all we do is redefine and reformulate familiar concepts and 
quantities as necessary. The earth and its environment are considered 
in their full complexity. Only at a much later stage do we intend to 
make approximations and only after a quantitative analysis of their 
effects. This paper which deals with the directional aspects of the 
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problem will be followed by subsequent ones which will treat the prob- 
lems of reference origins and scale, and also the question of how the 
reference frame can be established and maintained in practice. 

1. Fundamental Natural Vectors 

Natural vectors are defined as such by their property of being 
dependent only on some natural phenomena and consequently independent 
of any artifacts such as coordinate systems, reference models, etc. 
Consider a point P on the surface of the earth and another point Q which 
serves as a target being observed at some epoch T from the point P. For 
the epoch T we define a number of natural vectors at the point P, desig- 
nated as the fundamental vectors. 

Q - the Observational Vector . The light ray which travels from 
Q to P (or vice versa) is generally a space curve due to the refraction 
by the atmosphere. What we actually observe is the direction of the 
tangent to that space curve at the point P. This tangent line is 
defined as the observational fundamental vector and is denoted by Q. 

-r - the Local Vertical Vector . The gravity vector at the 
point P is denoted by T. Its magnitude is the value of gravity at P. 

We define the second fundamental vector - T, opposite in direction to T, 
to be referred to as the local vertical vector. 

- the Rotation Vector . Rotation is change of orientation of 
a body or mass element with respect to some inertial system. It can be 
found >y studying the space-like change of the velocity vector of mass 
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points with respect to inertial ?space. For example, if the space-like 
change is zero, that is, constant velocity at all points, there is no 
rotation, but only a translation. Let V be the velocity vector with 
respect to inertial space, then - rot V is by definition the rotation 
vector, also called the voruicity vector. Its magnitude is the instan- 
taneous rotation velocity. 

The definition separates reasonably rotation and deformation 
since the earth is not rigid, rof V just contains the antisymmetric 
part of the tensor grad V, whereas the symmetric part describes defor- 
mation. The earth rotation vector changes with respect to time due to 
precession, nutation and polar motion and with respect to space due to 
the deformability. 

X - the Ecliptic Normal . The ecliptic is the osculating plane of 
the space curve which the earth-moon barycenter is moving along. It is 
referred to a heliocentric system with inertial orientation. The vector 
X is the binormal vector of this curve. An approximation is the normal 
vector of the plane being spanned by the heliocenter (considered as 
fixed) and the earth-moon barycenter. 

Basic Angular Parameters 

Project the four fundamental vectors Q, - T, fi, and X onto a 
unit sphere centered at point P (Fig. 1). At any instant the four 
points are related by five basic angular parameters as follows: 

B altitude (observable) 

A azimuth (astronomically observable) 

$ latitude (astronomically observable) 
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H hour angle of vernal equinox 
E obliquity of the eclitic 



The vernal equine--* T is defined by ¥ = ft x X. The five angular para- 
meters depend on the positions of the four fundamental vectors. As the 
vectors were defined in general to be space and time variant, it fol- 
lows that the basic angular parameters are also space and time variant. 

2. Reference Model 

Analysis of a natural phenomenon is usually conducted through 
the introduction of an approximation, a so-called reference model. 

Using current knowledge of the phenomenon, a relatively simple model 
may be defined so that a reasonably good prediction of the phenomenon 
can be made for given space and time coordinates. In this section we 
define a reference model for the earth, the fundamental vectors, and the 
basic angular parameters defined in Section 1. 

Reference Model of the Earth . The model earth is defined dynam- 
ically (from the points of view of its gravity field and rotation) as a 
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rotationally symmetric level ellipsoid with major semiaxis a and eccen- 
tricity f- . The ellipsoid rotates versus inertial space with uniform 
velocity u> about an axis which is slightly inclined to its minor 
(figure) axis in accordance with a specified polar motion model. The 
mass of the ellipsoid m is equal to the mass of the earth, and the 
parameters a, e, and ui are selected so that the normal (model) gravity 
potential on its surface is constant and is equal to the gravity poten- 
tial on its surface is constant and is equal to the gravity potential 
on the geoid. The normal gravity potential at a given point, external 
to the ellipsoid, can be calculated from Gm, a, e, w, and the coordi- 
nates of the point where G is the Newtonian gravitational constant 
[Heiskanen and Moritz, 1967, pp. 64-67]. 

The orientation of the rotational axis versus inertial space 
for a given epoch is calculated by the currently adopted models and 
parameters of general precession and astronomic nutation. 

Geometrically, the model earth has a rigid irregular surface: 
the telluroid at a specified fundamental epoch [ibid., pp. 291-204]. Thus 
distances and angles between model surface points are assumed to be time 
invariant . 

The Fundamental Model Vectors . We define the fundamental, vec- 
tors of the model in a similar manner as for the natural case: 
q the observational vector is the straight line from the observing 
point P tc the target point Q as affected by aberration and 
parallax. 
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-Y the local vertical vector is opposite in direction to the vertical 
gradient of the normal gravity field at P 
to the model rotational vector at point P 
x the vector normal to the mean ecliptic plane 

The model fundamental vectors at a given epoch are related through »..<2el 
angular parameters similar to the natural ones as follows: 

3 model altitude 

a model azimuth 

<J> model latitude 

h model hour angle of the vernal equinox 

c obliquity of the mean ecliptic 

At a given epoch we can project on the unit sphere the natural 
and the model fundamental vectors as sha m in Fig. 2. The four dif- 
ferences 6q, fiy, fico, fix are called disturbance vectors. The disturb- 
ances in the basic angular parameters are 

fi(3 = B - 0 

fia - A - a 

5cf> = $ - <J> 

5h = H - h 

fie = E - E 

The mathematical relationships between the disturbance vectors and the 
disturbances in the angular parameters are given in Section 7. 


The fundamental vectors defined for the natural case and for the 
reference model vary in space and in time. The space-like variation of 
V is the difference between V + dV at a second point P + dP, in the 
neighborhood of P, and V at the same epoch. We can express the 
space-like variation of V as its partial derivative versus the space 
variable: 3V/3S. 

In a way similar to the space-like variation, we define the 
time-like variation of V at point P and epoch T as its partial deriva- 
tive versus the time variable: 3V/3T. The interpretation of the 

time-like variation is complicated by the necessity of defining the 
inertial frame as the common reference for the two states of the vector, 
i.e., and + ^ , . To simplify our treatment of variations, the 

fundamental vectors are placed in a hierarchy beginning with Q through 
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-r, n, X* up to 1 which is considered as any inertial vector. Thus, in 

order to obtain the absolute derivative of a fundamental vector V, ■=— , 

we have to differentiate both the coordinates of V and the base vectors 

defined by the fundamental vector on the next higher stage. Therefore ! 

we have to connect these base vectors with the inertial frame by means ! 

i 

of time-variable rotation matrices. These systems of base vectors and 
rotation matrices will be introduced in detail in the next chapter . 

In Table 1 we have listed certain phenomena causing the variations and j 

disturbances of the four fundamental vectors. A point to be kept in ! 


Funda- 

mental 

Vector 


Table 1 

Sources of Variations of the Fundamental Vectors 


Space-Like Variations 


Model 


Disturbances 


Time-Like Variations 


Model 


Disturbs .--s 


§ 


-r 


o 


parallax, 

aberration 


refraction 


relative motion 
of target 


perturbations in 
motion of target 


positional 

difference 


deflections of 
flie vertical 


local 

rotations 


constant spin 
rate, model 
polar motion 


correction to polar 
motion, spin rate 
variations, tides, 
mass redistributions 


luni-solar 
precession 
+ nutation 


correction to luni- 
solar precession 
+ nutation 


planetary 
pre cession 


correction to plane- 
tary precession, 
ecliptic wobble 


























mind is that certain phenomena associated with space-like variations of 
a vector are not necessarily time invariant and vice versa, as, for 
example, refraction or deflections of the vertical. 


4 . Natural and Model Triads 

The fundamental vectors defined in the preceding sections can 
Le used to define orthonormal vector bases, or triads. According to the 
vectors used there will be natural and model triads. 

Observational Triad - El . The three axes of the triad El at 

the point P and epoch T are defined by the vectors Q and -r as 

— 3 — 

El = norm Q 

El^ « norm [Q x (-I 1 )] 

— 1 — 3 — 2 

El - El x El 

Local Horizon Triad - E2 . The 
vectors -r and ft as follows : 

E2 3 = norm -F 

— 2 — 

E2 * norm [ft x 

— I — 3 — 9 

E2 = E2 x £2 

Equatorial Triad - E3 . The axes of E3 are defined by the vec- 
tors ft and X as follows : 

— 3 — 

E3 = norm ft 

E3^ ■> norm (ft x x) 

— 2 — 3 — 1 

E3 - E3 x E3 


axes of E2 are defined by the 


(-r)] 


li 


Ecliptic Triad - E4. The axes of E4 are defined by the vectors 


X and i. At this stage we introduce the inertial triad e which is a 
space- and time-invariant orthonormal vector base. Its specific orien- 
tation is not important at the moment and will be left undefined. Vec- 
tor i is parallel to axis e^ of the Inertial triad e. The definition of 
E4 is as follows: 

— 3 — 

E4 = norm X 

E4^ * norm (i x X) 

_ 2 — 3 — 1 

E4 = E4 x e4 

Note that axis 1 of E4 does not necessarily point towards the vernal 
equinox as is the case with the ecliptic system used in astronomy. 

The triads of the reference model are defined similarly, the 
only difference being the substitution of the model fundamental vectors 
q, -y, w, x for the natural ones. The model triads are denoted by 
lower case letters el, e2, etc.). 

The above definitions result in left-handed systems in El and 
E2, and in angular parameters (altitude, azimuth, etc.) in accordance 
with geodetic conventions (see [Mueller, 1969, pp. 32-42]). The triads, 
based on the same fundamental vectors, could also be defined more 
systematically (i.e., all right-handed), but in that case the angular 
parameters would not comply with presently accepted conventions. 

Transformation Between the Triads . We derive the orthogonal 
(rotational) transformations between the sequence of triads by intro- 
ducing three additional angular patameters, (see Fig. 3), 

which together with the basic angular parameters a, 3, <{>, h, and e 
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serve as parameters in the transformations. The sequence of transfor- 
mations is as follows: 

e4 - R 3 H' 3 ) e 

e3 ** R-^-e) R 3 <V> 1 ) e4 
e2 « p 1 r 2 (it/ 2 - 4>) R 3 (h)e3 
el - R 2 (tt/ 2-3) R 3 (a) e2 

where R^ (y) is a conventional rotational matrix around the j axis by 
an angle y (j = 1, 2, 3) [Mueller, 1969, pp. 43-44], 


is a permutation matrix of the axis k (k = 1, 2, 3) 



Fig. 3 


The transformations are orthogonal so the inverse relations can 
be obtained in general by reversing the order of the rotational matrices 
and also the sign of the rotational angle. The transformations between 
the sequence of the natural triads Ei (i = 1, 2, 3, 4) and e (the 


inertial triad) are the same, except that instead of the model angles 
one must use the natural parameters A, B, $, H, E and also 
^1’ ^2* ^3 ^‘ t ^ e ^ atter 8 rou P for fc he transformation between the inertial 
triad and E4) . 


5. Variations of a Triad 

Since the triads are defined by the fundamental vectors, it is 
obvious that their directional variations will involve a rotation of the 
triad. Such variations are possible in three dimensions: 1) in space, 

2) in time, 3) by the transition from the natural to the model funda- 
mental vectors or vice versa. 

Instead of analyzing separately the effects of these variations 

and disturbances, we shall study in a general way the influence of the 

variation of the fundamental vectors on the triads defined by them. It 

should be relatively easy, once the general formulae are available, to 

specify the kind of variation and the specific triad to which it 

applies. The same holds true for the disturbances . 

Let Z and D be two fundamental vectors (Z the "lower" and D the 

"upper" one). The triad of which D/ | D | ( = norm D) is the 3-vector is 

called E* = [E^ , E^ , E~* ]^. The triad of which norm Z is the 

— * 

3-vector is called E being defined as 

—3* — 

E = norm Z 


- 2 * 


= norm (D x z) 


— 1* —2* —3* 

E -EXE 


The representation of D in both systems is: 

— 1* 2* 3* — 

D = [D , D , ] E' 

1* 2- 3- . — | 

with the coordinates D = D =0,D = D , and, since the relation 
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between E* and E is 


E - R 2 (II/2 - <*>) R 3 (A) E* , 

— r 1* 2* 3*, -* 

D - [D , D , D ] E 


where [D 1 *, D 2 *, D 3 *] T - R 2 (n/2-$) R 3 <A) [D 1 *, D 2 *, D 3 ‘] T 

* [- [ D | cos 0, [ D I sin 4>] 


or 


D - -Id! cos $ E 1 * + I D ] sir. $ E 3 * . 


The fundamental vector Z is in the E -system: 


— 1* o* 3* _ 3 * 

Z - [Z , Z" , Z ] E 


■1 £ o 

with the coordinates Z = Z = 0, Z ■ Izl, 


or 


Z - i Z I E 3 * 


The variations of D and Z are 

dD = [dD 1 ' , dD 2 * , dD 3 ’] E* 

, 1* 2* 3* — * 

= [dD , dD , dD ] E 


with 

r i*n 

dD x 


- i . 

cos A sin $ dD 

2* 3» " 

+ sinAsin$dD - cos 4> dD 


2* 

du 

- 

- sin A dD 3 + 

- ,_2* 
cos A dD 


3* 

dD 


2* 

cos A cos $ dD 

2- 3* 

+ sinAcos$dD + sin $ dD 


dZ = [dZ 1 *, dZ 2 *, dZ 3 *] E 


Now let us construct the new base vectors E + dE after a small varia- 
tion of the fundamental vectors D and Z ^D + dD and Z + dZ. 
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—3* —3* — — 

E + dE ■ norm (Z + dZ) 


- norm {[dZ 1 *, dZ 2 *, i 3 * + dZ 3 *] E*} 

« rdZ dZ^ . g* 
l 37 » 3* , E 

Z Z 


— 2 * — 2 * 
E + dE 


norm [ (D + dD) x ( Z + dZ) ] 

norm { [dD 2 *(Z 3 * + dZ 3 *) - (D 3 * + dD 3 *)dZ 2 *] E 1 * + 

+ [(D 3 *+dD 3 *)dZ 1 *- (D 1 *^ 1 *) <Z 3 *+dZ 3 *) ]E 2 * + 
+ [(D 1 *+dD 1 *)dZ 2 *-dD 2 *dZ 1 *] E 3 *} 


[- 


dD 


2 * 


, 1 * 


3 * 2 * 2 * 

+ . 1 t 


1* 3* 
D x Z J 


E 1 * + dE 1- * - (E 2 * + dE 2 *) x (E 3 * + dE 3 *) 


[ 1 , 


dD 


2 * 


. 1 * 


3* 2* 

D dZ 

1* 3* * 
D Z 


1* 

“3*J E 


Collecting the new base vectors in one column matrix, we obtain 


1* — 1* 

E x + dE 

_ 2 * — 2 * 
E + dE 

—3* _3* 

E J + dE 


E + dE ■ (I + ft) E 


where the antisymmetric matrix 

r 


a 


2* 3* 2* 

dD D J dZ 

1* + 1* 3* 

D D Z 


dZ 


1* 

i3* 


dj£ P 3 *dZ 2 * 
D 1 *' D^Z 3 * 


dZ 


2 * 

3* 


dZ 


1* 


.3* 


dZ 


2 * 


, 3 * 


is the Cartan matrix ft [Grafarend, 1977, pp. 159-160]. Expressing the 
elements of the D-vector in terms of the E' -frame we get 


sinAsee$7=r 


dZ 


1* 


dD 2 ’ dZ 2 * 

- cosAsec<3>-p= — + tan$7 = 7-~ 

b Z 


ft 


-slnAsec*^ 


dD 


D 


dD 


+ cosAsec^^rssr; tan& 

D 


kdZ 


2 * 


dZ 


2 * 


dZ 


dZ‘ 


o* 


Now apply the general expressions derived above to any of the geodetic 
triads. For example, by Identification of D with ft and Z with -r , of 
dD with polar motion and dZ with a change of the vertical direction, we 
find the influence of polar motion and of a change of the vertical onto 
the orientation of the horizontal system. A and $ are then longitude 


and latitude, 



x and 



« -y the components of polar motion. 



and 


dZ 


2 * 


the angles of vertical change in north-south and 


east-west directions respectively. In order to get the horizontal sys- 
tem north-oriented, some signs have to be changed. Thus we finally 
obtain 
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— 1 — 2 — 3 

dE2 - (-sin A sec 4> x - cos A sec $ y + tan^k^) E2 + k, E2 

— 2 --1—3 

dE2 ■■ (sin A sec $ x + cos A sec $ y - tan 4> kj) E2 + k 2 E2 

dE2 3 » -k^ 1 - k 2 E2 2 

There are many similar applications of the general formula, for 

instance the influence of a change of the vertical and a motion of the 

target on the observational triad, or the dependence of the equat. 

system on planetary precession, luni-solar precession and nutation. 

While in the first two examples the motion is relative to an earth-fixed 

observer, it is described relative to an inertial system in the latter 

one. The general formula is valid for both cases. 

We can interpret the Cartan matrix as a rotation matrix of 

* 

three differential Cardan angles about the three E -axes, the first 

dZ 2 * dZ 1 * dD 1 ' 

angle being — — — , the second — — — and the third sin A sec 

1 z | | z ] |d| 

do 2 * dz 2 * i 

-cos A sec + tan $ — . Later on we shall call them T in the 

l°l M . 

first level (observational triad), V in the second level (horizontal 

1 i 

triad) , £ in the third level (equatorial triad) and y in the fourth 

level (ecliptic triad), i = 1,2,3. 

6. The Commutative Diagram of Triads 

To obtain a better insight into the interrelations between the 
various triads, we will construct three-dimensional structures to be 
referred to as the E(e) Towers or the Commutative Diagram of Triads 
(see Fig. 4). Each point in the diagram represents a certain triad 
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according to the ]abel attached to it. The straight lines between the 
points represent orthogonal (rotational) transformations between the 
respective triads. The overall organization of the diagram is as fol- 
lows : 


space-like 

variations 


time-like 

variations 


E-tower tower of the natural triads (solid lines) 

e-tower tower of the model triads (dashed lines) 

^evels 1, 2, 3, 4 - according to the type of triads, 

i.e., observational, horizontal, etc. 
the lines parallel to the space axis represent 
space-like variations of the triads 
the lines parallel to the time axis represent 
time-like variations of the triads 
disturbances the diagonal (dotted) lines which run on level i 
between an Ei triad and the corresponding ei 
triad. These are the only connections between 
the E-tower and the e-tower and represent the 
disturbances explained in Section 2. 

The diagram thus represents all triads , their space- and 
time .’.ike variations, and model disturbances at a single point P. In 
order to identify space- and time-like variations, we introduce two 
indices ( j , k) which follow the symbol Ei or ei of the triad. Both j 

and k can be 1 or 2, where index j - 1 stands for triads at P and j = 2 

at P -I- dP. In a similar manner k = 1 stands for triads at epoch T 

while k = 2 at T + dT. Thus the index (1, 1) indicates the situation at 

P at epoch T, (2, 1) is after a space-like, (1, 2) after a time-like 
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variation; (2, 2) represents the situation when both space- and 
time-like variations affected the trial (1, 1). 

Interlevel Transformations . In Section 4 we derived the inter- 
level transformations along a typical sequence Ei(j, k) , i * 1,2, 3, 4. 
Fig, 4 shows the pairs of parameters involved in a transformation 
between two adjacent triads along a column of the tower: A, B; H, $; 

etc. As these Interlevel parameters are space and time variant, it is 
obvious that they carry j and k indices matching the column of triads. 
We have identified the various parameters and the respective triad 
columns of the tower structure where they apply in Table 2. For com- 
pactness of representation, denote by a the vector of model angular 
parameters as follows : 

a T = [a, 3 , h, <p, e, ip 2 , 

Following the notation Introduced in Section 3, we denote space varia- 
tions by 30/9S, time variations by 9a/9T, disturbances (natural minus 
model) by <50, space variations of disturbances by 3(60)/3S, and time 
variations of disturbances by 9(<50)/9T. 
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ihr 






Table 2 


Interlevel Transformations 



Transformation, Parameters 

Tower Column 

Oi 

cr 

• 


ei a.i) 

c a 

3o* 

<r+ ^ ds 


ei (2,1) 

or a 

o 

+ t£ dT 

dT 


ei(l,2) 

O* 

, 3 o 3S 

tr+ SI dS + 3T dT 


ei (2,2) 

CTa 

to + cr 



Ei (1,1) 

Os 

to + o + 



Ei (2,1) 

07 

6 a + a 


dT 

Ei (1*2) 

Os 

to + o + 

'H + ^ 2 ] as+ [H + ^ 

dT 

Ei (2,2) 


Inlevel Transformations . We have defined in Section 5 the dif- 
ferential inlevel transformation vectors. In Fig. 5 we can see a total 
of 12 such vectors for level one. As the changes in the space and time 
variables dS and dT are differential and the diagram is commutative, 
there are seven independent conditions to be fulfilled: 



T 11 - -T 6 + T 10 + T 8 
*12 ' *7 + *11 * *9 
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space axis 


hr. r > 

and therefore only five independent t vectors loft. These represent 
the lol lowing varlat Iona In the triads of level ono: 

l j apaoo- 1 Iko j 

variat Iona of model triada 

i., time- like > 

T.j disturbances 
spare- like 

variations of natural triada 

t c time- like 

The various inlevel transformation parameters (t , v, etc.) ran he 
expressed aa functions of the relevant apace- and time- like variations 


of the fundamental vectors aa shown In Section 5, 




7 . Variations In the Basic Angular Parameters as a Function 
of Variations in the Fundamental Vactors 

We are faced with a large number of transformation parameters 
required to relate the various triads in the towers, we hcve already 
taken a step towards reducing their number in the inlevel transforma- 
tions. We will complete the reduction process by expressing the varia- 
tions of the basic angular parameters (interlevel transformations) as a 
function of variations of the fundamental vectors and show that the 
transformation between any two triads in the towers is de t mdent on the 
variations of the four fundamental vectors only. 

Following ideas in (Crafarend, 1977, pp. 207-212], in Fig. 6 
we have four triads which together form a closed loop of a commutative 
diagram. This loop is used as a typical example, and therefore the 
subscripts of t and v (which are v^» i.e., disturbances) are not 
indicated. From previous sections we have 


e2(l ,1) - R E (a,8) el (1 ,1) 

Kl(l.l) - R C (T) el(l,l) 

E2(l,l) - R (v) e2 (1 ,1) 

V» 

K2(l ,1) - R^a + 6a, 3 + 6P) F.l ( 1 , 1 ) 


where R^a.B) 


(— cx) ^(3 - tt/ 2), Eulerian rotation matrix 
R 3^ T 3^ R 2^ T 2^ R 1^ T 1^ and 


R^v) - R 3 (v 3 ) R 2 (v 9 ) R 1^ V 1^’ Cardanlan rotation matrices 

Tims 


R^V) e2(l ,1) - R^(a + 6a, p 4- 6?) El(l,l), or 

R c (v) R E (a,3) - ya + 6a, 3 + 53) R c (t) 
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E2 (I, I) 

* 

a+8a t 

£+8/? 


EKi, i) 


From here we can arrive, hv patient algebra, at the following expres- 
sions : 


flex 


-eosp 

0 

S i up 


0 

0 

-1 


- 




T + 






0 

-1 

0 


-s ina 

cosa 

0 


Treating in a similar way the other loops of the columns Ki(l,l), 
e i ( 1 , l ) , we get 


\Shl 


-C08t|' 0 

~sin<}> 


0 

0 

-1 " 

i 

m 



V + 




_dif> 


0 1 

o . 


-sin h 

cos h 

CL 


P'1-,1 


' 0 

-s inr 

I'OSf 


' 0 

0 

-1 

1 

m 




r. + 




. 5t: 



0 

0 


COSljij 

sini/'j 

0 


Substitute into r, \', f., p their equivalents, perform the multiplica- 
tions, and rearrange. The results are summarized in Table 3. The 
matrix presented is actually the matrix of partial derivatives of the 
basic angular parameters vs. variations of the fundamental vectors. It 
should be obvious that the matrix would not change if we considered 
space-like variations or time-like variations instead of the 




see 
































disturbances in the derivation as long as Fig. 4 is a commutative dia- 
gram. Table 3 represents the situation in the model. For the natural 
parameters, the relationships are, of course, identical. 
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Appendix A 


Differentials of a Compound Rotation Matrix 
Preliminaries 

Analytical expression for the differentials of an orthogonal 
matrix R which represents a sequence of elementary rotations is the 
subject of this Appendix. Rotation matrices are used in ortho- 

gonal coordinate transformations as shown in [Mueller, 1969, p. 43] 
where 6 is the angle of rotation and i is the axis about which the 
rotation is performed. 

The differentiation of a rotation matrix R^(6) with respect to 
the angle 0 is obtained by pre- or post-multiplying the R^(6) matrix by 
a skew symmetric matrix 

9R, (6) 

-55- - W e > ' 


The matrix is defined as the i layer of the skew-symmetric e.^^ sys* 
tern as shown in [Lucas, 1963]. The rotation and Lucas' matrices are 



"i 

0 

0 “ 


COS0 

0 

-sin0 


COS0 

sin0 

<f 

Rj_(e) - 

0 

COS0 

sin0 

R 2 0) * 

0 

1 

0 

r 3 (0)« 

-sin0 

COS0 

0 


_o 

-sin0 

COS0_ 


_sin0 

0 

cos0_ 


0 

0 

1_ 



"o 

0 

o" 


"o 

0 

-1 


'o 

1 

o“ 

L l“ 

0 

0 

1 

L 2 = 

0 

0 

0 

V 

-1 

0 

0 


_0 

-1 

0 _ 


_1 

0 

0 . 


_0 

0 
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A rotation matrix ^(9) or a product of two or more rotation mat- 
rices are orthogonal 3x3 matrices with the following two properties: 
(i) The determinant is equal to one. 

(ii) The inverse is equal to the transpose. 

The above properties of a 3x3 orthogonal matrix A can be utilized for 
deriving the elements of the adjoint matrix of A. As is well known the 
adjoint of a nonsingular matrix (the transposed matrix of its cofac- 
tors) divided by its determinant is equivalent to its inverse 



According to the properties of A as stated above, i.e., 
|A| * 1 and A ^ * A^ 
one has 

adj. A = A^ 
or explicitly 



(a 22 a 33 _a 32 a 23 ) 

“ (a 12 a 33 _a 13 a 32 ) (a i2 a 23 _a 22 a 13^ 


a ll 

a 21 

a 31 

adj A “ 

" <a 21 a 33“ a 3l a 23 ) 

(a li a 33~ a 31 a 13^ a ll a 23~ a l2 a i3 ) 

s 

a 12 

a 22 

S 32 


- (a 2i a 32“ a 22 a 31 ) 

~ (a ll a 32~ a l2 a 31 ) (a ll a 22~ a 12 a 21^ 


®13 

a 23 

a 33_ 


Use the above result in deriving an expression for the matrix product S 

S ** aba t 

where A is a 3 x 3 orthogonal matrix and B is a skew* symmetric matrix 
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®11 

fl 12 

a l3 


"0 

b 3 

1 

CM 

U3 

1 

A - 

a 21 

a 22 

a 23 

B - 

- b 3 

0 

b t 


- a 31 

a 32 

a 33- 


_ b 2 

- b l 

0 


Perform the multiplication, regroup to obtain 



(a 12 a 23 -a 22 a 13^ b l 

(a 

0 

+ (a 2l a 13“ a 23 a ll )b 2 

+ (a 


+ (a 22 a ll“ a 21 a 12 )b 3 

+ (a 


0 

(a 


33 12 32 13^ 1 

31 a 13 -a 33 a ll )b 2 
32 a ll“ a 31 a 12 )b 3 
33 a 22" a 32 a 23 )b l 
+ (a 31 a 23" a 33 a 21 )b 2 
+ ( a 3 2 a 2i“ a 3i a 22^ b 3 

skew-symmetric 0 

Using the property of the adjoint of an orthogonal matrix (A), 


0 (a 31 b l +a 32 b 2 +a 33 b 3 ) “ <a 21 b l +a 22 b 2 +a 23 b 3 ) 


(a ll b l +a 12 b 2 +a 13 b 3 ) 


skew-symmetric 


Differentials of a Sequence of Rotations 

The compound rotation matrix R which represents a sequence of ele- 
mentary rotations (0.) is defined as their product 

J J 


R = R. (e n ) ... R. (0JR, (0.) 
i n n i 2 2 i x 1 


Derive an expression for the partial derivative of R with respect to 
one of the angles 0^ where j = l,2,...,n and in a form which is con- 
venient for programming on a computer. 


30 


Partition R into three parts 


R - A R. (0, )B 

4 j 3 


where A, R, (0.) and B are orthogonal. 

j J 


L^B can be represented as BQ^ where is a skew symmetric matrix with 
elements which are a function of B. 


Q, - b\b 


Using the expression for S as developed earlier for each of the three 
cases i » 1,2,3 one gets 



■ ° b 13 

" b 12 


0 

b 23 

“ b 22 


0 b 33 

-»J2* 


- b l3 0 

b u 

Q 2* 

" b 23 

0 

b 21 

V 

-b 33 0 

b 31 


_ b 12 " b ll 

0 


_ b 22 

" b 21 

0_ 


_ b 32 " b 31 

0_ 


The resulting partial derivative of R is thus 


33J ' A * lj <V BQ i ' RQ i 


The variation of the R matrix as a function of variations of the 


0^ angles is obtained now easily from the above results 


SR 


3R 

30 


60 + 


n 


3R 

30 


j 


60. + 

j 


3R 

30. 


60. 


6R - R • Z Q ifl ■ R • i) 

j-1 j J 
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where Q Is a function of the i row of the product of the j - 1 ele- 
mentary rotation matrices to the right of R. (0.). 

J J 

Differentials of Cardanlan and Eulerlan Rotation Matrices 

There are two special types of compound rotation matrices which 
have been used extensively in deriving the various relationships in the 
E-tower : 


Cardanlan rotation matrix 
R c (a,0,Y) - R 3 (Y)R 2 (0)R 1 (o) 

and 

Eulerlan rotation matrix 
R E (a,$,Y) ■ R 3 (y)R 2 Ot/2-3)R 3 (oi) 


Using notation and formulae developed in the 
obtains for a Cardanlan matrix, 

" cosBcosY cososinY + slnasinBcosY 
R c « -cos0sinY cosacosY - sinasin0sinY 
_sin 3 - sinaco83 


preceding section one 

sinasinY - cosasin3cosY 
sinacosY + cosasinpsinY 
cosacos3 


8R C 3R C 3R 

15“ " R C L i ; W " R 3 (Y>R 2 <3)L 2 R 1 (a) ; jf - - L 3 R c 



0 0 0" 


~0 sina -cosa 


0 cosacosB sinacosB 

Q ■' 

001 

0 -1 0 . 
In* 

V 

-sina 0 C 

_co8a 0 0 . 

V 

-cosacos3 0 sin3 

_-slnacos3 -sin3 0 


6R c - R c • [ Va + Ve + Vy ] " Vc 
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“c 


0 sina6B + cosacosBSY -cosa6B + sinacos36y 
0 6a + sin36Y 

_skew symmetric 0 


The elements of the fl c matrix are differentially small, thus 
^ + 6^-^ • (I-4« c ) 

» Rq ( a , 3 , Y) Rg (6a+sin36Y , cosaS3-sinacos36Y , sina63+cosacos36Y) 

The derivation of a variational equation for the Eulerian matrix is as 
follows : 


*E 


" cosacosYsin3 - sinardnY 
-cosasinYsin3 - sinacosY 
_ cosacos3 


sinacosYsin3 + cosasinY 
-sinasinYsinB + cosacosY 
sinacosB 


-cosBcosY 

cosBsinY 

sinB 


8r e Sr e t 

ssr * y -3 ! - 5 T ■ v^> 






■ 0 

1 

0" 


' 0 

0 

cosef 


0 

sinB 

-sinacosB 

Q “ 
s a 

-1 

0 

0 

V 

0 

0 

sina 

V 

-sin3 

0 

cosacosB 


_o 

0 

0_ 


_-cosa 

-sina 

0 . 


sinacosB 

-cosacosB 

0 


l %f a + Vs + <W 


0 6a + sin36Y cosa63 - sinacos36Y 

0 sina53 + cosacos36Y 

Lskew symmetric 0 


33 


h +< w h< v 

- ^(OjBjY) *R c (»ina6g+co8aco866Y,-co»o60+»inoco8P6Y,6a+*in36Y) 


Note the similarities between the ftg and matrices: 
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Appendix B 


Differential Relationships Between Model and Natural 
Triads, Vectors and Angular Parameters 

Derivation of the differential relationships between model and 
natural quantities as presented in their final form in the main text are 
the subject of this Appendix. The results obtained in the last section 
of Appendix A are used extensively. For the sake of completeness, cer- 
tain formulae given in the main text are repeated. 


Levels 1 and 2 


el(l,l) Q 


a , $ 


e2(l,l) 6 



El(l,l) 


a+6a, 3 + 53 


E2 (1,1) 


Fig. 6.1 


The disturbances (6a, 63) of model azimuth and altitude respectively as 

Y 

well as the components of the two rotation vectors T » T ^ Tg]j 
T 

V = [v x \>2 Vj] are regarded as differentially small angles sp that the 
Cardanian rotation matrices R,(t) and R c (v) can be written as follows: 


R c (T) 
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From the commutative diagram in Fig. B-l, 


E2(l,l) -R 3 [ - (a+6a)]R 2 [(g+6g) - Tr/2JE1<1,1) 
E1(1,1)«R c (t) *1(1,1) 

E2(l,l) -R c (v) e2(l,l) 

e2(l,l) -R 3 (- a )/R 3 (e-Tr/2) *1(1,1) - R 12 *1(1,1) 


where R^ 2 


'singcosa -sina cosgcosa 
sings ina cosa cosgsina 
„ -cosg 0 sing . 


Using the formula for variation in R^ 2 as derived in Appendix A, 
E2(l,l) - R 12 (I+Q 12 )E1(1,1) 

where 


"o 

0 

-f 


0 -sing 

0" 


" 0 

-sing6a 

-6g “ 

0 

0 

0 

dig + 

sing 0 

cosg 

6a ■ 

singfia 

0 

cosg6a 

_1 

0 

0. 


0 -cosg 

0- 


. <$3 

-cosg6a 

0 _ 


From the four equations above and substituting the expressions for 6R^ 2 , 


r c (v).r 12 (M 12 )r c (t) r ; 2 


12 


t 3 - sing6a 


-t 3 + singda 

L T 2 + 6g 


-t 2 - 5g 
t 1 + cosgfia 


- cosg6a 


r; 


12 


R^ 2 is an orthogonal matrix and so the development for S in Appendix A 
can be applied: 
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" 1 v 3 - v 2 


-cosg(Xj+cos06a) -singsina(Xj+cosg6a)-cosa(x 2 +60) 



+sin0(x 3 -sin06a) -cosgsina(x 3 -sing6a) 

* v 3 1 V 1 

m 

1 singcosa(Xj+cos$6a)-sina(x 2 +60) 



+cos0cosa(x 3 -sin06a) 

_ v 2 " V 1 1 . 


skew symmetric 1 


from which it follows after regrouping: 


v. 


3J 


0 -sinct 
0 cosa 
-1 0 


' 



singcosa 

-sina 

cosgcosa 



6a 






• 

_6g_ 

4* 

sings ina 

cosa 

cosgsina 

• 

- 



. -cosg 

0 

sing 



L t 3 . 


The last expression in a compact notation is 


v 


12 


6a 

6b 


+ R 12 t 


Noting that * A.^ * 1 and also being orthogonal, the last 

T T 

expression premultiplied by and R^ 2 respectively yields: 


6a 1 
L6g 


A 12 v " A 12 R 12 t 


t * R 12 V - R i2 A 12 


f 6a 
63 J 


or explicitly 




r v ii 


r x ii 

'6a' 


o 

0 

1 

M 




"-cosg 0 sing" 




*s 


• 

V 2 

+ 


• 

x 2 

-6g- 


.-sina cosa 0 • 




. 0 -10. 
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" Y 


“-cosfS 0- 

'6a" 


-sinfJcosa sinpsina -cos0- 

T 2 

U.- 


0 -1 
- sin0 0- 


+ 

-sina cosa 0 

-cos0cosa cos|3sina sin$- 



Levels 2 and 3 

The same approach is followed in the derivation of differential 

expressions for levels 2-3 using the commutative diagram in Fig. B-2. 

The disturbances <5h, <5<f> of the respective model hour angle of vernal 

equinox and latitude are regarded as differentially small angles as 

T 

are the components of £ * £ 2 S-j]. The derivations are given with' 

out further comments. 



E2(l,l) 

h + 6h,<f> + 6<f> 
E3(l,l) 


E3(l,l) - P 1 R 3 (h+<Sh)R 2 (7r/2-(|)-6(J))E2(l,l) - R 2 3 (1+fi 23 )E2(1,1) 
E2(l,l) - R C (V) e2(l,l) 

E3(l,l) = R c (5) e3(l,l) 

e3(l,l) * P 1 R 3 (h)R 2 (Tr/2-<Ji) e2(l,l) - R 23 e2(l,l) 

where P^ is the permutation matrix for the reversal of the first axis 
[see Mueller, 1969, p. 43], 
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R 23 


sin<f>cos h -8 in h 

- sin<j>sin h cos h 
L coscf) 0 


C0S(f>C08 h 

cos<(>sin h 
sin<f> 


and 



- 0 

0 

1*1 


-0 

8in<t) 

0 - 


0 

sin<f>6h 

6 <p ~ 

fi 23 " 

0 

0 

0 

+ 

-e- 

<0 

-sin<() 

0 

COStf) 

6h « 

-sin<|>6h 

0 

coscf)6h 


■ -1 

0 

0- 


" 0 

-COS<f> 

0 - 


-- 6 <p 

-cos<{>6h 

0 ~ 


From the set of four equations above 


r 0 V 3 + sin<|> 6 h 


R c (£) - I + R 23 


0 


L skew symmetric 


+ &f> " 
+ cos<f>6h 



0 


Due to the permutation matrix the determinant of the 1^3 

T 

matrix is -1. Accordingly Ad j . - -R^. 

Skipping a few obvious steps the following is obtained: 


1 

ifn 


- 0 
0 

-sin 

cos 

h~ 

h 

•<Sh' 

+ 

- sin<J>cos h 
sin<{>sin h 

sin h 
-cos h 

-cos<()cos h- 
-cos<{>sin h 

■v f 

V 2 

k- 


L -i 

0 

- 



- -cos <|) 

0 

-sin(j) 

c 

l v 


5 


A 23 


■6h* 


23 


V 
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Premultiplying as in levels 1 and 2 and regrouping 


'6h' 


■ 0 0 

-r 

r ?i i 


• “COS<f> 

0 

-sin4>‘ 


-v 

- 64> J 

" 

--sin h cos h 

0- 

CM 

+ 

. 0 

1 

0 . 


V 2 





K 3~ 





Lv 3 J 


-v 

V 2 

m 

r- -COS <f> 0” 

0 1 

-6h- 

+ 

- sin<|>cos h sin<fisin h -cos<f>- 
sin h -cos h 0 

U,- 


-sin<p 0* 



— cos^cos h -cos<f>sin h -sin<f>- 



Levels 3 and 4 

As in the upper levels the disturbances 6\p^ and 6e are differen- 
tially small angles as are the components of the y rotation vector. 

The derivations are presented without comments. 



E3(l,l) 


+ e + 6e 


E4(l,l) 


E4(l,l) - R 3 (-^ 1 -6^ 1 )R 1 (e+6e)E3(l,l) - R 34 (I+fi 34 )E3(l,l) 
E4(l,l) = R c (y) e4(l,l) 

E3(l,l) - R c (C) e3(l,l) 

e4(l,l) = R 3 (-^ i )R 1 (e) e3(l,l) - R 34 e3(l,l) 

where 
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34 


p cos\p^ -£»inJ|)^cose -sinij^sine 

sinsf/^ cos^jcose cosi/^sine 
L 0 -sine cose ■ 


and 



O 

O 

0 

1 


- 0 -cose -sine 


- 0 -cose6i|^ -sineSty£ 

n 34' 

0 0 1 

6e + 

cose 0 0 

6^. 

cose&J^ 0 6e 


1 

0 

1 

t— * 

0 


-sine 0 0 - 


-sinefiiK -6e 0 ^ 


From the set of four equations above, 


V M > - i+v 


0 ^^-cosed^^ -^-sined^j-i 


h + 6e 


skew symmetric 
The resulting three matrix equations ares 
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"V 


r 0 

cosij;^- 

r«*,i 


-cosif^ 

-sini^cose 

-sirtj^sine- 


^2 

3C 

0 

sinij^ 

JL 

•6e ■ 

+ 

sinty^ 

cos^cose 

cosij^sine 




--1 

0 - 



- 0 

-sine 

cose 



V 



+ R 
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K 


r 6 M 


' 0 

0 

-1- 


~V 


• 0 

-sine 

cose' 


'*1 

1 

.<Se - 


- cossJj^ 

sini^ 

0. 

• 

^2 

+ 

.-1 

0 

0 . 

• 

=2 













- ?o- 
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-Kf 


- 0 

1-j 

^2 

- 

sine 

0 

-£,- 


--cose 

0- 


■6i|i 

.fie 


+ 


cosi{^ sinij^ 
-sintj^cose cos^cose 
--sin^sine cos^sine 


0 - 


-Vf 

sine 

• 

y 2 

cose- 


L v 


Combination and Summary of Differential Relationships 

In this section the formulae derived in the first three sections 
of Appendix B are combined with the results presented earlier. 

For compactness adopt the following notation: 



where q, ~Y, 
The matrix D 
to vector j . 
sions : 


w, x are the model fundamental vectors, 
stands for the partial derivative of vector i with respect 
Earlier we have derived the following differential expres- 


firi = D 1 * • M + D 71 • T 
V T 

6k = d| • £ + • v 

<$X = D X • u + D* • £ 

y 5 


V - D v • fin + D v • T 
n t 

£ * Dj* • 6k + • v 

y - dJJ • fix + r.^ • £ 


x = d t • fin + d t • v 
n v 

v - • fik + • £ 

£ = • 6x + ^ * y 

x y 


and the following: 
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T - Dj 6q + D* y 6(-Y) 

V " D ! y 6(-Y) + fiw 

£ ■ 5 w + fix 

(i) x 

Ji ■ fix 

x 

Substituting the second into the first group of equations 

fin - DVfiq + (dV y +dJ)_ v y >s(-y)+dV fiw 

5k - dV y 5(-y) + (D^+D^)fi U + fix 

<$X - D &) 5 fiu) + (D X D 5 +D X D y )fix 

^ 111 £ x y x 

Multiplication of the Dj matrices followed by rearrangement of terms 
results finally in Table 3 in the text. 


Appendix C 


Applications of the Differential Relationships 

Here examples of the application of some of the relationships 
presented in Table 3 are given. The examples have been selected from 
two important areas of analysis, i.e., space-like and time-like 
variations of the fundamental vectors and their effect on variations of 
the basic angular parameters. 

Space-like variations . Fig. C-l is the disturbance column of the 
E-tower over the first three levels. The disturbances, i.e., 
the differences between the natural and the model basic angular para- 
meters (a, 3, h, <J>) and fundamental vectors (q, -y, w) , are dif- 
ferentially small angles. The pairs of orthogonal components of 6q, 

6 (— Y) and have the following interpretations (see Table 1): 

<5q^ refraction in altitude 

Sq 2 refraction in azimuth 

6 (-y)^ meridional component of the deflection of the vertical 
6 (-y )2 prime vertical component of the deflection of the vertical 
6o)^ nonparallelity of the rotation axis in (ecliptic) longitude 
6(^2 nonparallelity of the rotation axis in obliquity 
From the first two rows of Table 3, 
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Thus the above equation relates possible errors (corrections) ' in refrac- 
tion, deflection of the vertical, and parallelity of the model (ellip- 
soidal) rotation axis versus natural rotation axis to those in azimuth 
and altitude. Assuming 6q to be zero and with a slightly different 
notation, we have the generalized Laplace conditions as shown in 
[Grafarend and Richter, 1977]. 

The above set could also be utilized directly as linearized obser- 
vational equations where 6a, 63 are the respective (observed minus 
model) azimuth and altitude and 6q, 6 (— y) and 6k) are the unknowns. 

Time-like variations . Fig. C-2 shows the time-like variational 
column of the F-tower at the second, third, and fourth levels. In this 
case we are considering natural angular parameters (H, $, E) and 
their values 6T later, denoted by H"*, and E". In accordance 

with Table 2, the expressions for the time-like variations of, e.g., 
the parameter H, are 

II' - H + || 6T 

- h + 6h + || 6T + 1^- 6T 

Subtract the corresponding model quantities and rearrange to obtain 
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(H' - H) - (h' - h) 




Similar notation is adopted for -y, w, and x. Now apply the third 

and fourth rows of Table 3 and write the following expressions: 

• • • 
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where 

6h variation in the disturbance (natural minus model) of the 
hour angle of vernal equinox 
6(f> variation in the disturbance of the latitude 
6(-y) effect of local (plate) motions plus differences (natural 
minus model) in polar motion and spin rate 
6to difference (natural minus model) in luni-solar precession 
and nutation 

Sx difference in planetary precession 

The above equation thus relates errors (corrections) in earth rotation 

(prr-cession, nutation, polar motion, spin rate) to those in latitude 

and hour angle (longitude + Greenwich sidereal time) . It could also be 

• • 

utilized as linearized observational equations where 6h, 6$ are the 
• • • 

observables and 6(-y), So), and <5x are the unknowns. 
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1. Introduction 


The fundamental vectors at a point on earth, their reference 
model, and the commutative tower of triads — the E tower — were intro- 
duced and studied in [Grafarend et al., 1979] as a first step in our 
investigation for reference frames in geodesy and geodynamics. Time 
was defined as the fourth independent coordinate, and an attempt was 
made to distinguish between natural — observable — quantities and their 
models corresponding to present day knowledge. The concepts presented 
there (vectors, triads, parameters, transformations, variations, etc.) 
were only directional. Distances, coordinates, linear velocities, scale 
and deformations were not considered. Consequently, no metric data of 
any kind could be analyzed with the help of the E-tower alone. 

In the following, we introduce the tower of origins (P-tower) 
which complements the directional E-tower in defining concepts, identi- 
fying parameters, and analyzing interrelationships and variations of 
positions and distances between points in space and time. The approach 
follows closely the one employed in [Grafarend et al., 1979]. The 
distances, coordinates, linear velocities of the various points, 
regarded as natural (real) quantities are paralleled by a set of models 
of the same in a one-to-one correspondence. As in the E-tower we are 
interested in the difference between the real and the model quantities 
to be represented subsequently as functions of a selected set of param- 
eters. The two towers are closely related in sharing certain concepts 
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and parameters and, in fact, it would have been impossible to present 
the P-tower without repeated reference to the E-tower. 

The tower of origins (P-tower) presented and studied in the fol- 
lowing chapters should not be regarded as a problem-solution-procedure 
type of report. It is rather an attempt to provide a method of analy- 
sis, to lay down foundations, to create a consistent and logical 
language and nomenclature for a subsequent study and solution of speci- 
fic problems. Many of the results in terms of concepts, relationships 
and variations may seem trivial and not necessarily new. But this is 
exactly the purpose of this report, i.e., to redefine, reorganize and 
systematize certain aspects of our present knowledge and understanding 
of the geometry, kinematics and dynamics of the earth without resorting 
to too many basic assumptions and hypotheses. We have tried to identify 
and clarify parameters and phenomena which apply to directions (E-tower) 
as well as to positions and distances (P-tower) between points in space 
and time. The creation of this common basis is essential for our 
future treatment of specific problems where we should be able to use as 
necessary a combination of concepts and formulae derived and associated 
with either of the two towers. 

The ultimate goal of our studies of reference frames for geodesy 
and geodynamics is the establishment of a conventional terrestrial 
coordinate system (CTCS) through the combined analysis of a selected set 
of high quality observations (laser ranging, radio interferometry, etc.). 
The CTCS should represent in some average but nonetheless well-defined 
manner the space-time behavior of the earth vs. inertial space. Dynamic 
or geometric variations of the earth in space and time would be referred 
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to the inertial frame of reference through the CTCS. The P-tower pre- 
sented in this report is another step toward the achievement of the 
above goal. 
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2. The Tower of Origins 


The relative positions and notions of points in space and time 
which are characterized as the origins of various reference frames for 
geodesy and geodynamics are presented and studied in a diagrammatic 
structure to be referred to as the tower of origins or the P-tower. 

The overall appearance, organization, and notation of the 
P-tower (see Fig. 1) are similar to those of the E-tower. The points 
in the diagram symbolize certain physically 'ingful points at a 
given epoch. Capital Pi denote natural-real origins, while their models 
are denoted by pi. The integer i signifies the level of the origin and 
assumes the values of 1, 2, 3, or 4, for the topocenter, bodycenter, 
barycenters, respectively. 

On a given level the points are organized along two axes: the 

space axis and the time axis. The integers within the parentheses 
(j , k) are the space and time indices of the point to be interpreted 
as follows: 

j - 1 point related to the observer 
j ■ 2 point related to the target 
k ■ 1 epoch T 

k * 2 "next” epoch T + dT where dT is a differentially small 
time interval. 

One should note the different interpretation given here to the j index 
as compared to the corresponding j index in the E-tower: In the 

P-tower Pi(l, k) and Fi(2, k) are two distinctly different (not adja- 
cent) points, which, in general, have different velocities in space. 

The level of a point depends on its nature and on its .unction which is 
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associated in general with the measurement of distances, directions or 
gravity. 

Topocentric (observational) level PI 

A point is considered at the topocentric level if it serves 
either as an observing point or as a target. Stars and quasars are not 
considered as target points since their three-dimensional coordinates 
are not known with equal precision. The principal point of a telescope, 
an EDM instrument, or of a radiotelescope are a few examples of 
observing points. The principal point of an artificial satellite's 
transponder or laser retro-reflector are a few examples of target 
points. 

Bodycrntric level P2 

The center of mass of a body serves as origin on the bodycentric 
level. A body is defined here as a conglomerate of mass points which 
are connected to each other fairly rigidly so that variations in rela- 
tive positions between the mass points (deformations) are small as com- 
pared to the overall size of the body. The earth and the moon are 
typical examples of such bodies and their respective centers of mass 
are points of the P2 level. We see that a P2 point has a definite 
physical meaning although it cannot be directly reached by observations. 
A point of the PI level is normally located on the surface of a body 
and as such is associated with a certain P2 point which is the same 

body's mass center. Exception to this rule is a close satellite of a 

planet (the earth or the moon) which is defined as a PI point while its 

P2 point is the mass renter of the planet. 
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Barycentric levels P3, P4 

A point is considered at the barycentric level If it is at the 
center of ms? of a set of bodies. The selection of the set is more or 
less arbitrary and thus identity of the P3 point depends on the composi- 
tion of the set (its elements) . There may be several barycentric levels 
according to some hierarchy. A good example for a P3 barycentric level 

(I) origin is the earth-moon barycenter. As the earth and the moon are 
a subset of the solar system set (the sun and the planets) we can define 
a P4 origin at the barycenter of the solar system (barycentric level 

(II) ). It should be obvious that one could continue with P5 at the 
barycenter of our galaxy, etc. 

Inertial level p 

The inertial origin p is defined as a point which is fixed or 
moving with uniform velocity in inertial space [see Goldstein, 1965], 

The positions and motions of all the points in the P-tower are referred 
to this p point in accordance with the lavs of Newtonian mechanics. 

The points in the diagram are marked either as full (black) 
circles or as hollow (white) circles depending on whether they repre- 
sent a natural point or its model. Thus, in Fig. 1 we can distinguish 
between the natural P-tower (the black points) and the model p-tower 
(the hollow circles). 

The lines between two points in the double tower represent vec- 
tors in inertial space. The interpretation of these vectors depends on 
the axis to which the vectors are parallel and also on the nature of 
the points being connected by it. 
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We will examine first vectors at the topocentric (observational) 
level. For example, let Pl(l, 1) be an observing point on the earth 

surface and PI (2, 1) represent a target on the lunar surface. As the k 

index (in the parentheses) is 1 for both points, the epoch T at which 
both points are defined is the same. 

The vector Pl(l, 1) PI (2, 1) (see Fig. 2) represents the natural 
geometric distance and direction between the two points. Analogously 
the vector Pl(l, 2) PI (2, 2) represents the natural distance and direc- 
tion between the same two points only at a "later" epoch T + dT. 

The vector which connects the positions of the same point at two 
different epochs (T and T + dT) is defined as the linear velocity vec- 
tor of that point vs. inertial space. For example, 

Pl(l, 1) Pl(l, 2) velocity of Pl(l, 1) at T 

Pl(2, 1) Pl(2, 2) velocity of PI (2, 1) at T 

The interpretation of the vectors connecting the points pl(l, 1), 
pl(2, 1), pl(l, 2), pl(2, 2) is the same as above but for the model. 

The differences between the instantaneous positions of the natural 
points and their models are represented by the following vectors (see 
Fig. 2): 

pl(l, 1) Pl(l, 1) = 3£l(l, 1) positional disturbance vector at 

epoch T for point Pl(l, 1) 

pl(2, 1) Pl(2, 1) = 6pl(2, 1) positional disturbance vector at 

epoch T for point PI (2, 1) 
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space axis 

Fig. 2 The topocentric level. 

pl(l, 2) Pl(l, 2) = 6pl(l, 2) positional disturbance vector 

at epoch T + dT for point Pl(l, 2) 
pl(2, 2) Pl(2, 2) = <Spl(2, 2) positional disturbance vector at 

epoch T + dT for point PI (2, 2) 
The vectors between the inertial point p and any of the natural 
or model points in the P-tower symbolize their position vectors in an 
inertial frame of reference with origin at p. By virtue of the above 
definition the P-tower is a cosanutative diagram of the vectors in 
inertial space, i.e., the sum of the vectors forming a closed loop is 
identically zero. Using the commutative property at the topocentric 
level (see Fig. 2) we derive the following relationships: 
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time axis 


P1(1,2)P1(2,2) - P1(1,1)P1(2,1) - P1(2,1)P1(2,2) - P1(1,1)P1(1,2) 


P1(1,1)P1(1,2) - pl(l,l)pl(l,2) - 6pl(l,2) - 6pl(l,l) 

An important property of a vector commutative diagram is that the vector 
relationships are independent of the coordinate system chosen to repre- 
sent those vectors. The components of the various vectors may change 
from one coordinate system to another; however, their magnitude as well 
as their relative orientation remains invariant. 

We will examine next a vertical wall of the P-tower (see Fig. 3). 
The k indices of all the points being 2 means that the wall represents 
a situation at epoch T + dT. In Fig. 3 we have used a shortened nota- 
tion for the vectors along the vertical lines as follows: 

PICT, 2) = P2(l,2) PI (1,2) 

P2(l,2) = P3(l,2) P2(l,2) 
etc . 

The interpretation of these vectors follows from the identity of the 
end points: 

PI (1,2) is geocentric (mass center) position vector of the 
observer at T + dT 

P2(l,2) is earth-moon barycentric position vector of the 
geocenter at T + dT 

etc . 

The vectors pi (1,2) and p2(l,2) are analogous to the above but for the 
model. 

The vectors which connect the model points with the corresponding 
natural points are defined as positional disturbances. For example, 
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Pl(l,2) - pi (1,2) - 6pl(l,2) - 6p2(l,2) 

tPlTl,2) + P2(T,2) ] - tpl(T,2) + P 2?I,2)] - 6iT?l,2) - 6p 3(1,2) 


etc. 

We will complete the examination of the structure and signifi- 
cance of the P tower by studying the interrelations between points on 
one of the time-variation walls as shown in Fig. 4. The meaning of the 

vectors connecting points along a column has been discussed above. The 
• » 

two vectors Pl(l,l) and P2(l,l) connecting the geocentric position 
vectors PI (1,1) and PI (1,2) are the respective linear velocity vectors 
vs. inertial space of the Pl(l,l) and P2(l,l) points at epoch T. Using 
the property of commutativity, we can write the following: 

PUT, 2) - P1(T,1) - Pltl,l) - Pltr.i) 

• • 

P2(l,2) - P2(l,l) - P2(l,l) - P3(l,l) 

etc. 

The expressions on the right-hand side represent the relative linear 
velocities of the observer vs. the geocenter and of the geocenter vs. 
the earth-moon barycenter, respectively. An interesting corollary is 
the following inequality: 

^ pttt.i) * pitr.D . 

Summarizing our discussion of the P-tower structure and the 
interpretation of the various points and vectors in it, we see that it 
can serve as a convenient means for representing and studying the whole 
range of positional and velocity information related to points in the 
natural world as well as in its model. 
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PI (1,1) 
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P3(l,2) 


P2(l,l) 


P3(l,2) 


P4(l,l) 


P4(l,l) 


-*#P4(1,2) 


Fig. 4 A vertical "time" wall of the P-tower 

It should be kept in mind that certain vectors in the P-tower can 
be null vectors due to the two end points being coincident. For 
example, if Pl(l,l) and Pl(2,l) are both points on the earth surface, 
the points P2(l,l) and P2(2,l) represent the same point, i.e., the geo- 
center, and therefore the vector P2(l,l) P2(2,l) is a null vector. For 
this case, we can easily deduce the following identities: 

P2d,l) = P2(2,l) 
p2U,l) = p2C2,l) 

6p2(l,l) = 6p2(2,l) 
etc. 
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3. Barycentric and Bodycentric Levels 


In Chapter 3 the nature and interrelationships o£ origins at the 
barycentric and bodycentric levels are studied. The major objective is 
to identify the positional disturbances and their time-like variations 
at these levels with inadequacies In current theories and respective 
constants. In particular we study the problem of possible dependence 
of second- and third-level disturbances on the rotation of the earth 
and its mass distribution. 

In Fig. 5 we have reproduced part of Level 3 of the P-tower 
relating it directly to the p point. Bypassing Level 4 in the above 
figure is the equivalent to the assumption that the solar system bary- 
center P4 and its model p4 are coincident and are taken as the inertial 
point. The vectors P3(l,l), P3(l,2) and their difference P3 (or, 

9 

equivalently, in this case P3) represent the motion of the earth-moon 
barycenter P3 with respect to the barycenter of the solar system. As 
the p3(l,l), p3(l,2), and p3 represent the model of the same, computable 
with current theory, it should be obvious that inadequacies in that 
theory will be represented by the respective disturbances. Accordingly, 
<5p3(l,l), 6p3(l,2) and their difference ^ 6p3 are all non-zero vectors. 

Little as we know at present about the 6p3 vector and its time- 
g 

like variation <$p3, we can at least state the following: The theory 

of motion of P3 about the barycenter of the solar system is a function 
of the combined masses of the earth and the moon, the masses of the sun 
and the other planets in addition to constants of integration (or zero 
epoch state vectors). Accordingly, phenomena such as (i) the motion of 
the mass centers of the earth and the moon vs. their barycenter F3, 
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P3 0,2) 



P3 earth-moon system barycenter 
p3 model earth-moon system barycenter 

p inertial point (identified with the solar system barycenter) 
Fig. 5 The Sarycentric level 

(ii) the mass distribution within the earth or the moon vs. their 
respective mass centers, (iii) the rotational motion in space of the 
earth or the moon, are not parts of the disturbances in the motion of 
the earth-moon barycenter. Another way of stating the above would be 
that measurements within the earth-moon system would not be sensitive to 
the <Sp3 disturbance or to its time-like variations. 

In Fig. 6 we have added Level 2 tc the previous case. F2 repre- 
sents the earth-moor, barycentric position vector of the geocenter (or 
selenocenter) F2. Using the commutative property of the loop formed by 

the four natural points, we can derive an expression for the vector 

* 

F2(1,1)P2(1,2) denoted in the diagram as P2 

• • 

P2 = P3 + P2 (1 ,2) - P2 (1,1) 
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Fig. 6 The bodycentric and barycentric levels 


but at Level 3 we had 

and so it follows 

3t 3t ■ 

Using the loop at level 2 and the above results, „e can writ 


3 «P2 = (P2 - p2) + A (pj . JJ, 


3T 


3T (<Sp2 “ 6 P 3 > - 4r (P2 - P 2) 



(P2 - p2) and its time derivative (P2 - p2) represent the difference 
between the real and the model motions of the geocenter about the 
earth-moon barycenter. We denote by c the ratio between the masses of 
the earth and the moon 


M 

C ■= rp ~ 81.3 
m 

and identify 

P2(l,l) as the position vector of the geocenter with respect to 
the earth moon barycenter P3 and 

P2(2,l) as the position vector of the selenocenter with respect 
to the same origin P3. 

From the definition of the barycenter of a two-body system we have the 
following (see Fig. 1): 

(i) P2(l,l), P2(2,l) and P2(1,1)P2(2,1) are collinear 

(ii) |P2(1,1)| + |P2(2,1)| - |P2(1,1)P2(2,1) | 

(iii) |P2(2,1)| / |P2(1,1)| - c. 


The above equations demonstrate the simple relationship between 
P2(l,l) and the lunar theory which in principle gives the components of 


P2(1,1)P2(2,1) and its time derivatives. 

. g ___ 

p2 and p2 can be computed from the current dynamic theory of 
oT 


earth-moon system (essentially the lunar theory) , while (P2 - p2) and 
its time derivatives represent corrections to that theory. 

From the above we can draw two conclusions which complement each 


other: 

(a) The positional disturbance of the geocenter 6p2 and its time 

g 

derivatives 7 ^ <Sp2 consist of two components which represent 
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corrections to the theories of motion of the barycenter of the 
earth-moon system and that of the geocenter (selenocenter) with 
respect to the barycenter of the solar system and to each other 
respectively. 

g 

(b) 6p2 and 6p2 do not depend on the rotation of the earth (or the 

moon) or on variations in its mass distribution. 

Summarizing this chapter and extending its conclusions to 6pl we 
can state: 

(a) Unaccounted perturbations in the theory of motion of the 
earth-moon barycenter with respect to the solar system barycenter 
dominate the 6p3 disturbances and their time derivatives. 

(b) Unaccounted perturbations in the lunar theory dominate the 
6p2 - 6p3 disturbances and their time derivatives. 

(c) The Level 1 disturbances 6pl or actually (6pl - 6p2) and their 
time derivatives are dominated by the rotation of the earth (or 
the moon) and by mass redistributions. 

It is important to realize that according to (b) above 6p2 does 
not have a diurnal motion and it is independent of inadequacies in the 
adopted gravity model of the earth (or the moon) . 
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4 . The Topocentric (Observer-Target) Level 


In this chapter we study problems which are related to the PI 
bodycentric position vector and its time-like variations. Our discus- 
sions are limited to Level 1 origins which are located on the surface 
of the earth (or the moon) . 

For a point on the earth (moon) surface there are three issues 
of fundamental importance which have to be carefully studied in order 
to understand the nature and significance of Level 1 positional dis- 
turbances and their time-like variations: 

(i) The rotational motion of PI with respect to an inertial frame 

of reference centered at P2. 

g 

(ii) The relationship between PI, ^ PI and the variable gravitational 
potential field of the earth (or the moon). 

(iii) The explicit definition (and realization) of pi — the model topo- 
centric origin. 

The proper order of introducing and studying the above three 
topics is not arbitrary as they are interdependent. Accordingly we 
begin with (i) proceed through (11) and finally end up with (iii). 

4.1. The Rotational Vector Its Model and Disturbances 

The time-like variations of the geocentric position vectors and 
their disturbances are strongly dependent on the rotational motion of 
the earth (moon) vs. its mass center. We vrill devote this sub-chapter 
to sharpening the concepts associated with the rotation vector 0 . for 
the earth and the rotational motion of PI around it. 

The rotational vector describes by its direction and magnitude 
the rotational motion of a point on the earth surface PI vs. the 
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geocenter P2 and with respect to inertial space. The time-like 
variation in position (linear velocity) of PI with respect to P2 is 
obtained by the well-known vector equation 


^ pi - n x pi 


which is rigorous for a rigid body. The motion of PI on the non-rigid 
earth vs. the geocenter P2 can be partitioned into two parts as fol- 
lows : 



n 


PI +[- 




where the first term represents variation in the direction of PI vs. 
inertial space and the second term is the variation in its magnitude. 
We will discuss in a subsequent sub-chapter the second term and its 
association with variations in the gravitational potential at PI. In 
the present sub-chapter we will be concerned only with the first term 
which is equivalent in form with the rigid rotation of PI about the 
mass center P2 (see Fig. 7). The absolute rotational motion of PI 
with respect to a non-rotating, inertial triad, represented in Fig. 7, 
by the inertial vector i, is quite complicated but can be partitioned 
into a sequence of simpler relative motions. For that purpose we 
define a sequence of rotational vectors - axes between i and PI ranked 
in the following order: 
i» ^g» PI 


where a higher rank is associated with the nearness to i (Fig. a) . The 
three rotational vectors are defined as follows: 
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ft g - is the spin vector. Its orientation vs. i is defined by the 

general (planetary + lunisolar) precession and by the forced terms 
of nutation. It does not contain terms of diurnal or higher fre- 
quencies. Its magnitude is changing in time with unpredictable 
variations to be determined by means of observations. 

ft_ - is the Eulerlan vector which rotates around f L with a nearly 

E S 

diurnal frequency and with a small angular amplitude |AQ |/|ft„| 

p o 

where Aft^ » ft^, - fl g . Both frequency and amplitude of Aft^ are 

unpredictable and can be determined only through observations. 

The A vector in magnitude and in orientation represents the 

polar motion phenomenon. The ft axis and its motion vs. i repre- 

E 

sents the complete solution of the differentia equations of 
rotational motion of the earth. 

Note: We should point out that both ft and ft are space invariant, 

b E 

i.e., they are the same in orientation and in magnitude for any 
point PI. Thus ft can be regarded as the instantaneous global 
rotational axis of the earth. 

ft - is the instantaneous rotational vector at PI. It has a nearly 

diurnal rotational motion around ft c and its angular distance from 


ftg is extremely small of the order of 10 seconds of arc. The 

Aft^ vector represents local motions of PI. 

In spite of the fact that ft or approximately ft^ are the true 

Instantaneous axes of rotation of PI, it has been demonstrated by 

Atkinson [1975] and also by Leick [1978] that the axis that can be 

detected directly by observations is the ft. axis, while the ft_, and the 

b E 

ft axes are unobservable. Intuitively, the above statement could be 
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explained by the fact that both and ft have a nearly diurnal rota- 
tional motion around ftg just like PI. An observer at PI cannot detect 

on a short time basis (one day or less) the motion of PI vs. ft_ and ft. 

E 

Only on a much longer time basis, i.e., days for ft_ and years for ft can 
one detect the accumulated effect of the small perturbations Aftp and 
Aft^ on the orientation of PI vs. ftg and through it vs. i. 

To recapitulate, we can state that the time-like variations of 
PI are defined by the following vector equation which is identical to 
the one written at the very beginning of this chapter. 



The above equation means, for example, that the angle between PI and 
ftg, the instantaneous geocentric colatitude, varies in time due to the 
combination of polar motion and local motions. In the second part of 
this sub-chapter we will define the model of the rotational vector ft. 

The dynamical model of the earth is defined as in [Grafarend et al. , 
1979] as a rotational level ellipsoid rotating with a constant spin rate. 
The orientation of its spin axis to with respect to an inertial frame 
is given by general precession and forced terms of nutation as pre- 
sently adopted. With respect to the axis of figure z (minor axis of the 
ellipsoid) the spin axis describes a cone with an amplitude 0'.'15 and a 
period of 1.1828 years. The sense of the model polar motion, thus 
defined, is counterclockwise as seen from the north. The two constants 
(O'.'IS, 1.1828 years) correspond to the average amplitude and period of 
polar motion between 1970 and 1976 [Markowitz, 1976]. 
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to 



Fig. 9 

The total rotational motion of the model geocentric position vec- 
tor pi is given by the following equation 

pl ■ (Wg + Aw p ) x pi « w x pi 

where Ato^ is the model polar motion vector. Its magnitude is 


A to 


2tt 0'.'15 _ 3 # 863085 10 " 6 rad/year, 


1.1828 p" 


it is normal to to , coplanar with to c and z (axis of figure of the ellip- 
S b 

sold) and points in a direction such that tog is between Aw^ and z (see 

Fig. 9). Ato rotates around to with an angular velocity slightly 
P s 

higher than | to | 


| tog J + ** 2306.4797 rad/year 


The angle between the Chandlerian axis to and the spin axis tog is 


-2- * 0'.'0003463 


to. 


74 



Fig. 10 


Because of polar motion in the model the instantaneous colatitude 
and the hour angle of the vernal equinox h of the position vector pi 
vary in time (see Fig. 10). These variations can be computed as fol- 
lows : 

0 » - | Aw | • sin (h - h ) 

P P 

h * |w | - | Aw | cos (h - h ) • cot 0 
S' 1 p 1 p 

As tht model of the earth is rigid the coordinates of pi in the x,y,z 
geocentric reference frame, fixed to the ellipsoid with z as the minor 
axis, are also invariant. After one polar motion cycle, i.e., 1.1828 
years, the 0 and h coordinates of pi will be back at their initial 
values. 
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In order to illustrate the feasibility of this type of parametri- 
zation of polar motion we evaluated the effect of the model polar motion 
on the geocentric equatorial coordinates of five stations. Over a 
period of 440 sidereal days which is slightly longer than the period of 
model polar motion (1.1828 years), we numerically integrated the time 
rates of the cr^ and h^ coordinates assuming spin and polar motion to be 
the only causes of their variation. The constants used were as fol- 
lows : 

to * 2301.1676 rad/year spin rate of the earth 
s 

hp m 2306.4797 rad/year spin rate of the polar motion vector. 

Atop * 3.863085 • 10 ^ rad/yaar polar motion magnitude 

The initial coordinates of the five stations are given in the 
following table: 

Table 1. Initial Equatorial Coordinates 


Station 

h 

a 

1 

10° 

50° 

2 

O 

00 

50° 

3 

154° 

o 

o 

4 

226° 

o 

O 

m 

5 

298° 

o 

o 

m 


The initial value of h was set at 180°. Table 2 shows the time-like 

P 

variations of the h, a coordinates of the five stations in arc seconds 
computed at 20 sidereal day intervals over a period of 440 days. 
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The values in the table are computed by subtracting from the numerically 
integrated h^, as affected by polar motion the equivalent h^, 
values without polar motion. 

In Fig. 11 we have plotted in addition to the varying a, h coordi- 
nates of the five stations also the varying position of a reference pole 
vs. the spin axis and the vernal equinox. The reference pole is defined 
in a way similar to that of the CIO pole, i.e., the angular distances to 
the five stations are invariant. 

We summarize this sub-chapter by writing up the equations for the 
disturbance vector in w, i.e., the difference between the real and the 
model instantaneous rotation vectors (at PI and pi respectively) as 
follows : 

610 » ft - w * (ft c -w c ) + (Aft -Aw ) + Aft„ * + 6w + < 5 w„ 

b S PP Jc S p £ 


where 

6u)g - is the disturbance in the spin vector of the earth, its first 

and second components being due to inadequacies of current theory 
and constants of precession and nutation and its third component 
representing spin rate variations. 

6w - is the polar motion disturbance vector; it is normal to and 
P S 

represents the di r £erence between real and model polar motions. 
< 5 ui£ = Aft^ - is the local component (space variant) of the 6w vector 
and is associated with local motions of PI. 


We repeat that 6u)_ and 6u) are global in nature, i.e., they are space 

bp 

invariant while is different, in general, for different points. 
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TABLE* '2 V* VARrATTONS~ITrEQU?nmiA'L*CnORDTNATES“DUr TO POLAR' MOTIGN 
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4.2 Gravity and the Position Vector 

The relationships between the geocentric position vector PI and 
its time-like variations on the one hand and the potential of gravity 
and its derivatives at PI on the other are studied in this sub-chapter. 

The time-like variations in PI were partitioned (see 4.1) into 
directional and magnitudinal components as follows: 

3T P 1 “ ^ x ?! + 3 t l P l) * |— i' 

01 jpi | 

We will study first the various causes for variations in the potential 
W at PI and their relationship to variations in the magnitude of the 
geocentric position vector. 

The gravity potential of the earth at PI (which is a point on the 
earth burface) is evaluated by the following well-known formula (see 
Fig. 12): 

W - G / £ dv + i (Q x PI) • (fixH) 

V x i 

where dv is an element of volume, 

p,H are the density of dv and its distance from PI respectively, 

G is the gravitational constant. 

The integration is extended over V which includes in our case the solid 

earth, the oceans and the atmosphere. W thus obtained would be the 

measured value of W from which the potential, of extraterrestrial masses, 

the tidal potential has been subtracted. Considering the total mass 

within V to be invariant 

/ pdv = const ■ m 
V 
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Fig. 12 



0 


we can differentiate W with respect to time to obtain the time-like 
variations in W 


9W 

3T 


G / 
V 



1 

£ 


U . 

9T 



where denotes variations in rotational potential. 

The first term in the integral is not associated with any changes in 
the relative distances between PI and other material points including 
possible target points P1(2,K) or other observing points P1(1,K). In 
other words, time-like variation in gravitational potential due to 
density redistribution within the earth are not accompanied by time-like 
variations in relative position. However, the position of the mass 
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center F2 does change (mass center shift) within the surface S and con- 
sequently the geocentric vector PI also changes. The variations in the 
position vector PI due to mass center shift are space invariant, i.e., 

they are the same for all the points on the earth. 

3W 

The second term in •%=; is explicitly associated with variations in 


relative distances between PI and the totality of mass points which con- 
stitute the earth. The phenomenon which dominates this term is the 
local differential motion of PI. One out of several causes for local 
motions is the elastic response of the earth to variations in the tidal 
potential. 

By definition, the "horizontal" component of motion of PI is nor- 
mal to -T (the direction of the local vertical) and produces zero vari- 

8P1 

ation in the potential. Thus, the only component of which is 
3W 

related to ^=7 is the vertical component, i.e., 


3P1 

3T 





where -r is the unit vector along the local vertical (see Fig. 13) and 
che approximation is permissible due to the small angle between PI and 


-r. 


The time-like variations in magnitude of the geocentric position 
vector PI are related to variations in the potential W by a modifica- 
tion of the well-known formula [Heiskanen and Moritz, 1967] which 
relates potential and height differences: 


3W 

3T 


- -g 
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Fig. 13 


from which it follows directly 


_ 3 _ 

9T 




3W 

9T 


In the second part of this sub-chapter we study the relationship 
between the local gravity vector -r at PI and the geocentric position 
vector PI, their models and the corresponding disturbances. 

Fig. 14 shows a schematic spatial diagram of the geocentric posi- 
tion vector PI; the -T local vertical vector ar.d their respective models 
pi and -y. The disturbances 6p2 and 6pl as well as 6(-y) are also shown, 
z is the axis of figure of the reference ellipsoid. We derive first an 
expression for the angle between -y and pi. 
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Fig. 14 


A right-handed Cartesian coordinate system x, y, z is defined 

(see Fig. 15) which has its origin at p2 and which is fixed to the 

reference ellipsoid. The ellipsoid, or equivalently, the x,y,z system, 

rotates vs. inertial space around the spin axis w which is inclined by 

S 

0V15 vs. z. The point pi is defined in the x,y,z system by its three 
Cartesian coordinates pl(x,y,z) or by the three geocentric spherical 
coordinates p, a', X which are the geocentric radial distance, colati- 
tude and longitude respectively. Since the model is assumed to be 
rigid he three coordinates are constant. 



z 


X 


(0 


S 




Fig. 15 


We will use also ellipsoidal coordinates of pi in the x,y,z 
coordinate system, namely u, 3» A which are convenient in that the 
gravity (normal) potential U of the model ellipsoid and its derivatives 
can be represented in closed formulae. 

For the computation of U and the components of its gradient y we 
will assume that the component to of to_ along z is equal in magnitude 

b 

to |to c |. The difference between to and |<o c | divided by to is negligible 
b * b 

-12 

- of the order of 10 . The value of the normal potential U at the 

point pi is computed as a function of the four parameters of the level 

ellipsoid a, e, m, to and the components of the position vector pi. 

According to [Heiskanen and Moritz, 1967] and utilizing ellipsoidal 

coordinates u, 3, A, the following expressions hold: 
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_ /p 2 - a 2 e 2 / 22 2 

“ -J 22 — 2 + V 8 e P 


„ .0 C08 CT 

$ “ arc sin *■ 

u 


0000 2 2 \ 2 

2 2_2 2_ t , (p - a e ) 

1 e p cos o + •*£-> ^ 


the Inverse relationship being 

„ \l 2 x J2 27 

p « \u + a e cos 8 


p ' * arc tan ^1 + ^ “ jf cot 8 

Longitude is the same in spherical or ellipsoidal coordinates. The 
potential at pl(u,$,A) is computed by the following formula (see ibid.) 
where the effect of the small equatorial (x-y plane) component of Wg 
has been neglected: 

TT , a\ G ® _ ae .1 2 2 q , . 2 0 1. , 1 2, 2 X 2 2. 2, 

U(u,8) * — arc tan — + •=- u) a (sin 8 - t) w (u +a e )cos 3 
ae u 2 qo 3 2 


where 


I [( 1 + 3 - !t) 

u a e 


ae 

arc tan 3 

u ae 


“1 
le J 


qo 


■i[( 


(1 - e j 

2 

e 


) 


1 + 3 ) arc tan — — — ■ - 3 

Vl-e 2 


Vl-e 2 


The vector along the gradient of U at pi is y the normal gravity vector. 
Its components along the ellipsoidal coordinates are: 

9U _ ~Gm 


9u , 2 2 2, 

(u +a e ) 


. ,,2 2 q , 1 w 2 a 2 / . 2„ 1\ 

+ w u cos 8 + — — — / sin 8 - jj 


qo 


3u ae 
o - arc tan — 
z z u 


a e 


3u 2 + 2a 2 e 2 
ae(u 2 + a 2 e 2 ) J 
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■ 57 f “ W^sin3 cos3 [a^ - (u^ + a^e^)] 

dp qo 


3U n 
8X " 0 

In order to obtain the components of y in the x,y,z system we need the 
Jacobian of transformation from u, 3, X into x,y,z coordinates. 

The independence of U from X Implies that the three vectors z, pi and y 
are coplanar. Accordingly, instead of transforming from u, 3> X into 
x,y,z, we transform from u, 3 into r, z where r, z, X are the cylindri- 


cal coordinates of pi and 
The transformation equations are simple 

J 2 2 2 Q 

r = \u + a e cos3 


V 2 2 

x + y is the distance from the z axis. 


z ■ u sin3 

The components of y in the r,z system would be computed then in a row 
vector form 


r 8U 8U q r 8U - 

L 87 87-* LaTT aftJ J 


8u 83- 


where 


J = 


!(I1 \lu 2 + 

HI) 


2 2 2 . 
u + a e sin 3 


u cos 3 


- sin3 


\l u^ + a t e‘ sin3 


2 2 
i e 

u cos3 

J 2 . 2 2 

Vu + a e 


The angle between y and the z axis is thus given by 


0 * = arc tan 


r 9U -« 

8r 
SU 

L 87 J 


87 


while the corresponding angle of the pi vector is (see Fig. 16): 

cot0 


* arc tan — * arc tan 


\fx? + a 2 e 2 



We point out that $ the complement of is different from the conven- 
tional geodetic latitude which is computed for a point on the ellip- 
soidal surface. We shall see in a subsequent sub-chapter that the 

the ellipsoid. The angle between -y and pi is 

and a few kilometers height above the ellipsoid 

10 ’. 


model of PI is not on 
thus 

Aa’ « 0* - a' 

Y pl 

and for mid latitudes 
it is of the order of 
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We will derive now an approximate expression for the angle between 
PI and the -r vector as a function of Act', the positional [<Sp2, 6pl) 
and the angular [<S(-y)] disturbances. The angle between PI and pi 
denoted by A is evaluated from the dot product: 

« cos A 

|P1| Ip1| 

From Fig. 14 we have 

PI * pi + (Spl - <Sp2) 
cosA * pi • pi + (?pl - 

|pi| | P i| 

a - r | pi | + (<Spl - 6p2) • 

|pi| L |pl| 

We expand cosA (A is a small angle) and regroup 

4 = - [ML + . -5-1 

\ L jpi| |pi| |pi| J 

If (<SpT - <5p2) is collinear with pi, A will be zero. A will be 
maximum for |pl| = |pT| from which we derive finally (see Fig. 17) 

~ l?l| 

If we define pi and p2 so that the magnitude of the difference 
(6pl-6p2)is of the order of a few kilometers, A will be a small angle 
of the order of tens of seconds of arc. 

Considering now that (-D , (-y) , PI and pi are not necessarily 
coplanar we can see according to Fig. 18 .hat the angle between PI 
and -P can be approximated by Acr’ the error being smaller than the sum 
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(6pl - 6 P 2 ) PI 


Pi 




\ 

\ 



Fig. 17 



[A + 6 (— y) J - The difference between the real AE' and its model Ac', 
i.e., the disturbance in Act' depends on the deflections of the vertical 
<5(-y) and on the positional disturbance difference (6pl - 6p2) . 


4.3. Time-Like Variations in Level 1 Positional Disturbances 

In Chapter 3 of this report we studied the general nature of the 
second- and third-level positional disturbances and their variations with 
time. In the case of Level 1 positional disturbances and their 
time-like variations we will be more specific. In this sub-chapter we 
derive the differential equations of Level 1 positional disturbances in 
terms of disturbances of the rotational vector ft and also in terus of 
variations in the magnitude of the geocentric vector PI. Assume that 
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the real geocenter P2 and its model p2 are coincident, i.e., the second 
level positional disturbance 6p2 is identically a zero vector. 

Later we will relax this condition and will show the resulting 
implications. The time-like variations of the geocentric vector PI, 

3pT * 

or PI was partitioned above (see 4.1, 4.2) into a variation in 
direction and a variation in magnitude (see Fig. 19). For completeness 
we rederive the expression for PI in a slightly different form: 



•— ■ — - • I PI I * Q x PI ■ (to x 6to) x (pi + 6pl) * 

|5r| _ _ _ 

to x pT+ to x 6pT+ 6to x pi + <5w x 6pT 



Fig. 19 
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The time-like variation of the pi vector by definition consists 
of a directional component only as the magnitude of pi is invariant 
(The model of the earth is assumed to be rigid.) 

3-r- — —r- 

pi ■ 0) X pi 

where u is the sum of the diurnal rotational (spin) vector w and Aw 

s p 

the model polar motion rotation vector. 

The expression for the Level 1 positional disturbance simplified 
by the assumption dp2 * 0 is 
dpi * PI - pi 

The time-like variation of dpi is obtained by differentiation as 
follows: 


ST 6pl " ST P1 " ST pl 


We substitute expressions derived above, neglect two terms of the second 
order (Atop x 6pl) , (dto x dpi) , and obtain 


S jTT 

6pl ' Tiii 


ST 


PI + dto x pl + to x 6pl 


Rearranging and substituting for Pl its equivalent we obtain the final 


form of dpi 


ST 


dpi * to X dpi + 6td X pl + -Jr I pl + dpi) (p*- + - - ^ - P — 

S (7 1 I _ I r> i I 


I pl + dpi I 


This is a set of three first-order differential equations of the posi- 
tional disturbance dpi with dto as an independent parameter. It demon- 
strates the relationship between the Level 1 positional disturbances and 
those of the rotational vector. If the positional disturbance dpi at 
some initial epoch is known, we can integrate numerically the differen- 
tial equations of dpi using rotational vector disturbances and 
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variations in the nagnitude of PI (or equivalently variations in the 
potential W) , whi^h have been determined from observations. 

We will show now that the above equation holds also for the case 
where <5p2 is not zero. Using a portion of the P tower as in Fig. 20 we 
can use the commutative property to derive the following (Level 4 is 
excluded without loss of generality) : 


3 

-gf 6p3 


6p3(l,2) - 6p3(l,l) “ H 1 - I? 1 


3T 


<5p2 - 6p2(l,2) - 6p2(l,l) 


9 .-=■ , 9P2 3p2 

3T 6p3 + 9T" ' 9T 


9T 


6pT « 6pT(l,2) - 6pT(l,l) “ gY 6p5" + fi x PI + 


3IP11 

9T 


PI 

IpTI 


- u x pi = 6p? + o) x PI + 6d) x PI - u x pi 


. il ri 

9T 


PI 

IpTI 


Neglecting second -order terma, substituting W g for a) and 


regrouping we have: 


•^•(6pl - <5p2) - W g x (6pl - 6p2) 4- 6w x pi 

0 s Ip 1 ! 


.-^|pl+(6pl - 6p2) 


The resulting vector differential equation of the geocentric positional 
disturbance (6pT - 6p2) is similar to the one obtained earlier for 

6p2 = 0. 

In the second part of this sub-chapter we will develop a specific 
model of the Level 1 origins for the earth. The models of all the PI 
points of the. earth constitute thus the geometrical model of the earth 
surface. ^ u rotational level ellipsoid discussed earlier in this 
report is the dynamical model of the earth. By making the distinction 
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between geometrical and dynamical models of the earth we actually define pi 
as a point which is not located on the surface of the ellipsoid . As we shall 
see, pi is located on the telluroid as defined by Hirvonen [1960] and as 
described also -in [Heiskanen and V'-ritz, 1967], 

The selection of the telluroid as the geometrical model of the 
earth is essential for establishing a clear and unambiguous relationship 
between the potential of model (normal) gravity at pi and its deriva- 
tives on the one hand and the geocentric model position vector pi and 
model gravity vector y at pi on the other. 

In the P tower we have denoted a point on the earth surface as 
either Pl(l,k) or Pl(2,k) where indices 1 or 2 indicate an observing 
station or a target at the epoch (k = 1,2). In the following dis- 
cussions we will drop the indices for convenience, as it will become 
clear that the particular values of the two indices within brackets are 
irrelevant. The fundamental vectors, discussed in [Graf arend et al., 1979], 
-r, are specifically referred to the PI point, where in particular F 
is the direction of gravity at PI and is parallel to the axis around 
which the geocentric vector PI rotates with respect to inertial space 
(see sub-chapter 4.1). 

We denoted the model of PI as pi and denoted the vector dif- 
ference plPl as Spl the positional disturbance of pi. Just as for PI 
above, the models of the fundamental vectors -y, to refer to the pi 
point. In particular the y vector is defined as the direction of model 
(normal) gravity at pi and u) is parallel to the axis around which the 
model geocentric vector pi = p2pl rotates in inertial space. In order 
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to focus on pi, PI and 6pl alone we will make the assumption that p2 
and P2 coincide, i.e., 6p2 * 0 (see Fig. 21). 

In principle <5pl cannot be and remain a zero vector due to the 
essential difference in the rotational motion of the two vectors pi and 
PI. ft and u are different in direction and in magnitude; |pl| is con- 
stant by definition (rigidity) while |Pl| varies in time due to various 
causes like tides, mass redistributions, regional uplifts, etc. 

We will define now the relationship between PI and pi (see 
Fig. 22) through the concepts c’> the height anomaly ? and the telluroid 
as described in [Helskanen and Mo; xtz, 1967] . In addition to the basic 
angular parameters ®, H which define the orientation of -r versus ft 
we introduce the gravity potential W at PI or actually the potential 
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Fig. 22 

difference W DO - W between the geoid and PI, The parameters of the 
reference ellipsoid (a, e, w, m) are chosen so that the model (normal) 
potential on its surface U 0Q is equal to W OQ . In Fig. 22 the quanti- 
ties 0, H and W define the position of PI in space and the direction 
of -r ther; . Apply to -T the disturbance S(-y) with an opposite sign 
to obtain (except for a small correction Ac") the -y vector. 

Beginning from PI we measure the height anomaly t, along the -y 
vector and obtain the pi point, i.e., the model of the PI point. The 
positional disturbance vector Spl thus is defined in magnitude by Z, and 
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by -y in direction. The point pi so defined is on the telluroid. 
According to the definitions of € and the telluroid the normal gravity 
potential at pi is equal to the natural potential V at PI or equiva- 
lently the respective potential differences versus the ellipsoid and 
the geoid are equal 


U 

oo 



W 


oo 



Now apply to ft the disturbance 6u) with opposite sign and obtain the 
direction of u) in space. Using a rigorous transformation from -y to 


pi (see 4.2) obtain the direction of pi in space. 


The magnitude of pi is obtained from -J at pi, the a angle between 
pi and z (o' after being corrected for model polar motion) and the para- 
meters of the ellipsoid. Thus, we arrive finally at the p2 point, the 
mass center of the ellipsoid. 

We can summarize in concept the above relationships as follows: 

(i) Three quantities (x, y, z) or (<f> , X, U) , are needed to define 
the pi and -y vectors. 

(ii) Two disturbances (<S(-y) and £) are needed to transform from pi 


and -y into PI and -I 1 . The two disturbances are represented by 
three numbers: two for 6(-y), the deflections of the vertical, 

and one for C, the height anomaly which is close in value to the 
undulations of the geoid. 

The disturbances <S(-y) and C as defined above correspond to the 
quantities which would be evaluated through well-known techniques of 
physical geodesy [Heiskanen and Moritz, 1967]. 

There are two basic difficulties involved in the above definition 
of the positional disturbance 6pl: 
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(i) The gravity potential W at PI varies in time and correspondingly 
the model potential U at pi which is equal to W should vary also. 
This, however, would require a non-rigid telluroid which contra- 
dicts our definition of a rigid model of the earth. 

(ii) and to are different in magnitude and in direction. As the 
above vectors represent the rotational velocity vectors of PI 
and pi respectively, it is obvious that the two points will not 
remain aligned along the -y vector, except at an initial epoch. 

A possible solution which allows us to retain some of the obvious 
advantages of the telluroid as the geometrical model of the earth with- 
out sacrificing the rigidity principle is as follows: 

The geometrical model of the earth is assumed to be rigid. It is 
defined as the telluroid at a specified zero epoch. From the zero epoch 
and on the positional disturbances vary according to the differential 
equation derived in the first part of this sub-chapter. 

4.4 Time-Like Variations of the Distance 
Between Two Earth Surface Points 

In this sub-chapter we study variations in the distance between 
points at the topocentric level in order to identify the global and 
local parameters which can be recovered. Consider the distance between 
two points on the earth surface, i.e., Pl(l,l) and Pl(2,l), the 
observing and the target points at the topocentric level of the P tower 
(see Fig. 23). As both points are defined on the earth surface their 
body-centric reference point P2 is the same (the geocenter) for both and 
so the vectors P2(1,1)P2(2,1) , P2(1,2)P2(2,2) , etc. are all null vec- 
tors. 
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We will simplify the notation in this sub-chapter and adopt 
the following: 

C - Pl(2,l) - Pl(l,l) 

c - iT(2,l) - iT(l,l) 

6c » 6pl(2,l) - 6pl(l,l) 
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where C, c, 6c are the respective observer-target vector, its 
model and its disturbance. From the commutative properties of the 
P tower (topocentric level) we can easily derive the following (see 
Fig. 23): 

C * c + 6c 

The rates of change of the 6bo”e vectors are reflected in the dif- 
ferences C(2) - C(l), c(2) - c(l) 6c(2) - 6c(l) and can be obtained 
by formal differentiation vs. the time variable : 

• * • 

C • c + 6c 

By C we will denote the rate of change of the magnitude (length) of the 
vector C. From sub-chapter 4.1 we have 

3 ~r — 

*§j£ pi - u x pi 

which when applied to the difference pl(2,l) - pl(l,l), and remembering 

that to is space invariant, results in: 

• 

C * 0) x c 

Note that c « 0, i.e., the distance between any two model points is 
Invariant according to the assumption of rigidity. 

The disturbance in the rotation vector 6oo is presented in two compo- 
nents as follows (see sub-chapter 4.1): 


6co *= 6a)- 4- 6oo 

g S p 

global component 

St 

local component. 


From sub-chapter 4.3 we have 

(6pl - 6p2) *= w x (6pl - 6p2) + 6u x PI + t~|P1|* 
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which when applied to the difference 6pl(2,l) - 6pl(l,l) results in: 

6 c a u x 6c + 6w g x C + [6(0^2 x pT(2,l) - 65^ x pl(l,l)J 

+ |pT( 2,1) I . - 3 |PT(1,1) | • ] 

T |pi(2 »i) j 9T |pT(l,l) I 

<5c“ = CO x 6^c + 5u XC+C6L + 6M) 

s 

T ?ith the above the rate of change of |c|, i.e,, C is 

• • • 

r s C • C m (c + 6c) • (c + 6c) 

C C 

* ^ [ii) x c • c+coxc • 6c + co x 6c • c + (J x 6c • 6c 
+ 6u x C • C + (6L + 6M) • C] 

o 

The first, fourth and fifth terms in the square brackets are zero 
due to the fact that two of the three vectors in the mixed vector pro- 
duct are the same. The second and the third terms cancel being of the 
same magnitude and opposite sign. Thus finally, we have the following: 

C = I • (6L + 6M) « - • (6L + 6M) 

C c 

Explicitly written the result is 


JL 

3T 


|pi(2,i) - Pl(l,l) | « 

|pl(2,l) - pl(l,l) | 
*| [ 6^£ 2 x pT(2,l)-6co^ 1 x pT(l,l) j 

+ ^ |fr<2 ’ i)l • J8£r & |E<1 ' 1> 


pid»i) 

Ipia.i) 


From inspection of the above equation we can state the following: 
(a) The rate of change of the distance between two earth surface 
points is independent of global phenomena. 
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Pl(2, 



PI (1,1) 


(b) The vector sum (SL + 6M) represents the difference in loc al 
horizontal and vertical motions (the relative motion) between 
the two points. 

In the last part of this sub-chapter we will study the effect of 
a shift of the geocenter (due to mass redistributions) on the distance 
between t\co surface points. 

Denote the shift of P2 vs. Pl(l,l) and Pl(2,l) by A and decompose 
it into three vector components A^, A^ along the directions of 
Pl(l,l), Pl(2,l) and Pl(2,l) x PI (1,1) respectively (see Fig. 24). 
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The component is normal to the plane defined by PI (1,1) and 
PI (2,1) and also to the vector C. Accordingly its contribution to the 
sum(6L + 6M) is also normal to C and so the dot product is zero: 

C • (SL + <5M) a “ 0 

The effect of A^ on the sum (6 l + 6M) can be represented by the equiva- 
lent parallel shifts of Pl(l,l) and Pl(2,l) in the opposite direction. 
The magnitudes of 6L and 6M due to A^ are as follows (A^ * |"S^ | ) : 

|6l| ■ A^ sin^ - 0 

^|PT(1,1)| - A x 

■^■|P1(2,1)| - A^ cos ip 

|6m| -Va/ ~ A^ cos^iJj * A^ sinifi 

The magnitudes of <5L and 6M being the same and by Inspection of Fig. 23 
we get finally 

6L + 6M * 0 

A similar proof can be derived for 

Thus we see that although the shift of the geocenter, A causes 
local variations in the orientation and the magnitude of PI vectors, it 
has no effect on the distance between PI (surface) points. 
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INVESTIGATIONS ON THE HIERARCHY OF REFERENCE FRAMES 
IN GEODESY AND GEODYNAMICS 


PART III: SCALE SYSTEMS: THE S'-TOWER 


by 

Erik W. Grafarend 
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The third hierarchic structure in Euclid space; 
the tower of geodetic scale systems 


0. Introduction 


While the hierarchic structures which rule orientation 
and origin (rotational and translational degrees of freedom) 
have been presented with respect to space-time geodesy in 
E. Grafarend (1978 a, b) and E. Grafarend, I. Mueller, H. Papo 
and B. Richter (1979), the third hierarchic structure wixl be 
introduced here, namely scale. Any vector space is furnished 
with the topological notion of length, here the lengths 
of geodetic reference vectors like the length of the gravity 
vector, of the rotation vector, of the ecliptic normal, etc. 
Beside directional parallelism scale parallelism is needed, a 
notion introduced by H. Wcyl (1952 p. 121-138). 

Spacelike and timelike changes of fundamental geodetic length 
with respect to a fixed length or scale unit (unit length, unit 
time and others) will be studied, extending first results of 
refraction studies in E. Grafarend (1976) where Weyl-geometry 
was used. The variations will be finally applied to the three 
base vector system (r, fi, T) which establishes three-dimensional 
geodesy. As a special technique polar and singular value decompo- 
sition are used in order to separate angular and dilatational 
distortions. The results can be embedded into the general theory 
of deformations introduced by C. Boucher (1978). 


r>r» 

« l* 


• if *4 * a 


FAGS EL.AMX MOT 
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1. The local structure of the scale system 


Prom the differential point of view two derivations of the 
basic scale structure in Euclid space are given. The relation 
to IVeyl geometry is emphasized. 

1.1 

Here, let us introduce v(x,y,z,t), a four-dimensional or 

space-time vector field which is a function of space-time 
12 3 4 

coordinates x -x, x =y, x =z, x =t in Euclid space. Hie vector y; 
is represented twofold, firstly with respect to an orthonormal 
triad (c^e^e^) such that its coordinates are (0,0, v) where 
v is the length of the vector, secondly, with respect to an 
orthonormal triad (e^o ,ey- >e-o) which is fixed in space-time 
or invariant with respeot to a translation in space-time. The 
base vectors are related by a rotation, e -> e = Re , where R 
is a threedimensional rotation matrix. Space - and/or timelike 
variations are studied by differentiation: 


1(1) v = (0,0, v) 

/V 



(0,0, sv 0 ) R 


j3° 



e i 

e 1 

Q 


1 1 

’ <NJ tO 

( -8W, 

1 (2) dv = (0,0, dv) 

fyj 

+ (0,0, v) 


t/ *7 

^3 


r _ 

1 


r~ 


d{ (0,0, sv Q )} 


R 


£1° 

^2° 

®3° 


+ (0,0,sv o ) d{R 


5l° 

S 2 0 

®3° 


} = 


{0,0,dsv Q + s dv Q )} R 


— — 

e 1 o 

e 1 

|2° 

+ (0,0, sv Q ) 

{dR 

e 1 0 

e 1 

~2° 

+ R 

do. 

a ^ 

A~ 

G70 
^ f) 



~ 3 ° 


^3 
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The length of the vector v has been expressed by the 
product of a a.'u.V s and a /aadwitu.' !r>ujth v . 

For exajnple, a length of 10 m is the product of the scale 
factor s = 10 and the fundamental length v * 1 m. In addition 
to the translational invariance of the /.>>< iv/.uv.w 

OjVtf e 0 we will assume that the fundamental length v 
is invariant with inspect to translation, too. Thus beside 
the postulate of directional parallelism d£^ = Q we have the 
of a. •>;/.’ /.'//it’tj dv = 0. These postulates lead 

to a variation of the vector field v given by 


1(3) dv = (0,0, dv) 



► (0,0, v) 



i (0,0, dsv) ) R 


( (0,0, dss" v) 


°l 

e , 

c: 


+ lO.O.sv lidR 


el 

5? j 


el 

e“ 


r e 

} ^ I 


+ (0,0, v) dRR* 


ol 

e“ 


or 

1(4) dv = ds s“'v 
1(5) do = dRlf 1 e 

Note that R is cm orthogonal matrix, |R| = +1, or R~ * = R' . 

A verbal fomulation of the fundamental result is this: The 
length of a vector v is changed under directional and scale 
parallelism proportional to the change of scale factor and 
the length itself, but inverse proportional to the scale 
factor. The orientation of the reference system e ( is changed 
under directional and scale parallelism proportional to the 
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change of directional parameters within the rotation 
matrix R and the base vectors e y itself, hut inverse pro 
portional to the rotation matrix. 

Fig. 1 illustrates the degrees of freedom of type trans 
lation, rotation, scale or origin, orientation, scale. 


Parallel transport of directions 
iC|c, e , Q1 e.o)' and length unit v i 



p(x,y.z.M p(x*dx,y*dy,z»dz,H 

or 

p(x,y,z,t*dH 



1.2 


Another derivation of the fundamental differential equations 
1(4), 1(5) originates directly from the group of transform- 
ations. According to Fig. 2 let us denote by v(x o ,y o ,z o ,t o ) 
a vector at a space-time point x,y,z,t. Both vectors coincide 
if we change orientation and scale by 

1 ( 6 ) v(x,y,z,t) = s R v(x 0 ,y o » z 0 »V 

1 (7) dv = dv e, 3 + v de^ = (ds s ^ + dR R by 
or 

dv = ds s - ^ v 
de = dR R ^ e 

F-ig. 5 illustrates the different postulates of parallel 
transport of directions and scale. 


Fi. 7 . 


° f fr f edom of translation, 

rotation, scale ’ 
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Fuj, Directions and scale under translational 
invariance: di^ « 0, dv Q « 0 


translation 
> 


scale at 

parallel 

point 

transported 

P(*o.y 0 . 2 0 . r o> 

scale at 


point 


plx.y.z.t) 

e.g. 

« 9 

1m 

1 m 



directions at 


point 

p(x 0 ,yo,Zo.t 0 ) 


translation 

> 



parallel 
transported 
directions at 
point 

p(x,y,z,tl 


1.3 

Hie classical treatment of length variation in differential 

2 2 

geometry is based on the quadratic fonn v = \ |v| | of the 
vector v. dv and dv are obviously related by 

1 (8) dv 2 = 2v dv 
1 ( 9 ) dv = -2-dv 2 
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leading to 


U10) 


T d ||v ||“ 


V. 


d In s 


II y li 


91ns 

9X 1 * * 


dx 1 . 


2. Hie global structure of scale systems 

Prom the integral point of vieu a derivation of the basic 
scale structure in Euclid space is given. The invariance of 
observables under the group of transformations is emphasized. 

2. 1 


Here let us introduce two vectors v(x o ,y o , z Q , t Q ) and 
,v.(x,y,z,t) at space-time points * 0 , y o ,2 o ,1: o and x,y,z,t, 
respectively, which are parallel under a translation. Both 
vectors coincide if we change orientation and scale by 


2(1) 

v(x,y,z,t) = s R v(x o ,y o ,z o ,t o ) 


2(2) 

« s y(x,y,z,t) = v(x o +6x,y 0 +fiy,z o +6z,t o ) 

- V .(WV 

2(3) 

fi t y(x f y,2,t) = v(x o ,y o ,z o ,t o +6t) 

- xtWV 


6 g v; is called spaaelike variation, { v, timelike variation. 

Let us introduce the rotation parameters by 


2(4) R = ^,*,0) = R 3 (0) R 2 (t -*) R 3 (A) 
2(5) R e (A+6A,$+S$,0) = R e (A,$,0) 


1 

-<$A 

-cos A6$> 


+6A 

1 

-sin A<5$ 


+cos A <5$ 
+sin A6$ 

1 _ 


+ 




2 


where '-‘ 7 indicates terms of second order. 


115 


2(6) 6v * 6s s ^ V + Rg(A,$,0) 6A Rg(A,4>,0) V 
where the antisymmetric matrix - A can be represented by 


2(7) 6A - 


0 

-6A 

-COS A64> 


+6 A +COS A6<J 

0 +sin A5 <t> 

-sin A6$ 0 _ 


2(8) n = RgfA.^.O) 6A Rg(A, 4>,0) » 


n 

-6A sin $ 
-6$ 


+6A sin $ 
0 

+6 A COS $ 


+6$ 

-6A COS $ 
0 


2(9) 6v » 6s s V 
2(10) 6e - ne 


We will prove next that positional angles and lengths ratios 
are invariant with respect to the underlying similarity 
transformation 


2(11) v **• Tv 

|V r ' 


s R v + t 

♦ »* r * 


2 ( 12 ) 


<T(y 2 -y 1 ) , T(y 3 -v 1 )> 

!l T (y. 2 -Xi)|lll T (y3- v .i)ll 


s 2 (y7-y-[) 

s2 llVXi ] IIX3-X1 II 


<v 2 - v 1> v 3~ v 1 > 

: v 2“ v 1i ; i! v 3" v 1 q - e - d 
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A 




H 

>7 

X 
\ — / 

s /(v,-y,)' R'R (v^v^' 

x2-n 1 

^1 T(v 3~v i | 

s /(v^-y^' R’R 

| v ^- v :i 


Fig. •/ is an illustration of the invariance of positional 
angles and lengths ratios in a space-time triangle. Related 
commutative diagrams for translation, rotation and scale are 
given in Fig. b 


3. 1-xainples 

Threedimensional geodesy will he based on three base vectors, 
namely [f , n, Tj ' , located at the topocentre and referring to 
the vector fields of gravity, rotation and eliptic normal. The 
base vectors are neither orthogonal nor normalized. The gravity 
vector determines the local A'ertical. The rotation field is 
constructed from the inertial velocity vector v of the topocentre 
by vorticity SI = rot v changing in space and time due to plate 
rotations and the dynamics of the planetary system. The eliptic 
normal is defined by the binormal vector of the curve of the mass 
centre of the earth in inertial space. The base vectors will be 
referred to a base vector system at initial epoch zero and space 
point zero, in detail by 


r 


mm 

To 


Vo 



*s/ 


r* 

n 

= RU 

n» 

A> 

II 

% 

SI® 

% 






which corresponds to a systematic set-up of type 


3(2) 


v = RU v = VR v . 

ivO 
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Fig. 4: Space- time triangle 



(B(x, y,z, t+6t) 
or 

e(x+6x,y+6y, z+6 z,H 


v(x,y,z,t+6f-) 

or 

v(x+6x,y+6y,z+6z,tt 


p(x,y,z,U6n 

or 

p(x+6x,y+6y,z+6z, t) 


Fig. 5: Commutative diagrams for degrees of 

freedom of type rotation, scale and 
translation 



It includes the polar decomposition (Cauchy decomposition) 
where R is a rotation matrix (|R| * + 1), U and V are right 
a>ui left stretch matrices being syranetric. The matrices are 
related by 


3(3) V = RUR'«=* RYR - U. 

A singular value decomposition of the stretch matrices is 

3(4) v = r v v*r; 

3(5) U = R U*R' 

v u u 

where 

3(e) V* * diag (Vj , v,, v~) 

3(3) U* = diag (Uj , u,, u-) 

and Vj, v,, and Uj, u,, u^ are eigenvalues. 

3(8) v = RR U*R' v » R V*R'R v 

-v U U V V ''O 

leads to variations of type spacelike and/or timelike 

3(9) dv = (dRU + RdU)y o = (dVR + VdR)y Q = 

(dRR' + RdlJ U" 1 R')v = (dW _1 + VdRR'V _1 )v 

or 

3(10) dv = (dRR' + R(dR IMR* + R U*dR' + R dU*R')R U* -1 R'R}v 

V U U U U U U U U 

= { (dR v V^r + RyVW + R v dV^)R v V^ -1 R^ + VdRR’V -1 }v 
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3(11) dU » diyj'R/ + lyJ'dR^ + R^U'R^ 

3(12) dV - dRyV*^ + R v V*dR^ r + lyFVRJ 

3(13) dv « dR R' v + (RcUyi'R^ + RR^JW) R u diag(l , t; , ~ )R’R^ 

12 3 

+ RR U diag(du r du 2 , du 3 ) diag (1 , ^ JR^R'^ 

"1 2 3 

3(14) dc * (dR v diagfv^, v 2 , Vj)RJ + \ diag (v 1 , v 2 , v 3 )dR^,}R v 

JiaR 4,’ ? 2 > v, 1R v i' + R v dia 8 (dv 1' dv Z- dv 3’ 

dia « 4,’?,4,^ tvdRR ' v ''i 

The tensors 

3(15) C = U 2 , B = V 2 

will he calledjng7;/- and left deformation matrices (Cauchy-Green 
matrices) which can be represented by 

3(16) C = r u U * 2r u = R u diag Cu 1 2 ’ u 2 2 ’ u 3 2)R u 
3(17) B * \V*\ = R v diag (v^, v., 2 , v 3 2 )R; 

What is the sense of all these strange computations? 

At first we have rotated the three base vactors by a proper 
rotation matrix R. Secondly we have stretched the three base 
vectors by the matrices U and V, respectively. The singular 
value decomposition allows the separation of angular and scale 
distortion. By the matrices R u and R y , respectively, we have 
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rotated the matrices U and V, respectively, into their 
principal directions. Along the principal directions there 
is only a change in scale of the three base vectors 
(J, g, ' . Thus we have found a decomposition into shear 
and dilatation , the off- and diagonal elements of the 
deformation matrices if we use this terminology. In general, 
the space-time change of geodetic base vectors can therefore 
be understood as a change in origin (translation), orientation 
(rotation) and scale. Fixed or translational invariant is al- 
ways the base vector system (j 0 , g Q , Y ]. In geodetic applicat- 
ions, the nine elements which describe the space-time change 
of a triplet of base vectors is parameterized in a slightly 
different way: The base vector g of rotation is projected onto 
the plane rectangular to the base vector T; the direction is 
called south. Orthogonal to south within this plane we direct 
east, equivalently by the vector product fiA-^T; the normalized 
triad as the final product is called the horizontal one. By a 
similar process applied to land we arrive at the equatorial 
triad. Angular parameters which connect these triads are always 
of type longitude and latitude. Totally there are six angles which 
connect the system of base vectors £r, fi, 4j] ' , which span the 
geodetic three-dimensional Euclid space locally, and the one 
[r , g Q , Tj . In addition, there is a space-time change of lengths 
||r|,|| gll, ||T|| parameterized by three saale factors referring to 
a fixed length system jlrjl , IlgJ , (lyji. 
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